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ABSTRACT 


This analysis presents a theoretical boundary-layer formulation 
which is valid for steady, two-dimensional and axisymmetric free-convection 
flow of any Newtonian fluid along surfaces having smooth, continuous temper- 
ature distributions. In addition, the stability of these flows is analysed 
theoretically using the small disturbance theory. The linearized forms 
of the disturbance equations are presented for two forms of small distur- 
bances: two-dimensional disturbances of the form of Tollmien-Schlichting 
waves, and three-dimensional disturbances of the form of Taylor-Goertler 
roll vortices. 

Numerical solutions are presented for steady, constant-—property, 
boundary-layer free-convection flow of air along long, inclined, isothermal, 
plane surfaces. These steady- or mean-flow solutions are subsequently used 
in the stability analysis to obtain numerical solutions for points on the 
neutral stability curves for both forms of the disturbances. 

The mean-flow solutions are found to be in good agreement with 
previous theoretical and experimental results, where these results are 
comparable. It is concluded from the mean-flow results that over the range 
of surface inclinations from 0° to 120° the most significant change in the 
profiles occurs within the range from 0° to 60° with only very small changes 
over the range from 60° to 120°. 

The results of the stability analysis are in good agreement with 
previous theoretical results, but not with previous experimental results. 
However, the difference between the theoretical and experimental results 
is attributed to an amplification process. The stability of the flow is 


strongly influenced by the surface inclination over the range from 30° to 


 poktelemre? reywl-qrebnnod teotsexoadd , adaseor eset ater 


paeicaeag sue’ obey tomorye £28 | bits Lanokefiontb~ows Ndiaiae: jane seca ay 


ui 


~reqnes reovatinoo <dtcome ynived aanmtrua: goole binl? neknosyoi we to 
beeylene et ewolt sass to yatitdsse oft) Cae teties ne -edoisudiases} & 


am2o% bestrsenbt sdT Kross Son aye Disme 543 yakew ie 


Ln 


—xutekb Iieme Yo aniot owd TOL betases1g S26 aia iin sonsdisebb § Sie 
| os ahineiaeaanmicinedtel to wrod ada to ssonsdauterb Lsnolenamtb=ows aga 
je. 

1913x909~r01 ys to. aro? os ad apomnday seth tude tenaelh bdaanene! bos 4 

: 


.egotsiov 


Loe 
e 
sift 


a i 
q 
FT 


eYIISgoTg~Itstanos .vbests rot alginate: exe endtiuloe Lata 
Asmaarizoe! ,bentiont ,amol yrois ats to wolt no Fiasvabo=aen seyel-Vae 
beeu vlanoupsadue S16 anokiuloe woli-nzem 10 ~ybsste sesnT veogatran) 

ori no aintog ot enotjuloe Iso lzoimun angen ow aleybens were ree odd 


.eeonsdxutakb odd Yo emrot dtod 10% Bevin vititdars texsia 


1 " 


daiw jnemeptgs boog mi od o2 brvot 818 enokauios woL nam odT ty. 


e7s etivesx sas) si1edw ved Lees iesnomtsagxo bas fsotséxosils ed 


egme7 93 x9vo aad) aptuter woLt~soon gis nose BebeSonnaial 4 ta 


ptt ok sgnado danodttngte team ada OSS ot °O mori , anokjantioat die kaert 
ee 
togneds ‘ema YI9v dao ria "0 03. °O) mx? sgos7 ed? ‘bist e193. #9 


bt y <n 


0% 
sce 


ry 
ew 


A 


r , ; . -08t ‘oar “08 mor , SgaBt pwr 


oho ye ‘cn. of i a ={I : : 4 + 
9 ee ems & ekectans aime? ‘atts og ealuest sit 


5 


- 
Oe oul 
De | iA 


AA : 7) va 


seaiuens Fatesernageé: meres iatw A Rassias baie aes ed | 


wi 
a mre a “” aati 
, . wa LIST i “| e 


sei i ; 
, 3 ; 


— 


ib 


120° with the most significant change occurring at about 85°. For angles 


below 85°, the results reveal that the Taylor-Goertler roll-vortex distur- 
bances determine the stability of the flow, while above 85°, the Tollmien- 


Schlichting wave disturbances determine the stability. 
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normalized radius of axisymmetric body 
normalized coefficient of thermal conductivity 
coefficient of thermal conductivity 


a reference length in the X direction representing the 
maximum distance from the leading edge at which the flow 
is examined. 


normalized coefficient of viscosity 
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of bulk viscosity 
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transformed coordinate 
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transformed pressure in [,& coordinate system 
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Subscripts 


c characteristic value; conjugate when used with vector 
notation, i.e. OV 


d departure 

- imaginary part of a complex number 

fo) hydrostatic 
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Lg constant pressure 

c real part of a complex number 

w wall 

(eines partial differentiation with respect to Z,n,& 
oo ambient, far from surface 


Superscripts 


associated with bulk viscosity 
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CHAPTER I 


INTRODUCTION 


1.1 Previous Analyses of Steady, Laminar, Free-Convection Boundary- 
Layer Flows. 


One of the earliest analyses of free-convection flows dates 
back to the theoretical work of Lorenz ile in 1881. He considered the 
problem of a free-convection flow along a heated vertical plate, and in 
his analysis, Lorenz assumed that the streamlines and the isotherms are 
parallel to the surface of the plate. However, it is now known that 
these assumptions do not agree with experimental observations. 

In 1930, Schmidt and Beckmann [2] suggested that the boundary- 
layer assumptions might be applied to the equations for at least some 
free-convection flows. Using these assumptions, they obtained solutions 
for a free-convection flow of air over a heated vertical plate. Schmidt 
and Beckmann also performed an experimental analysis for the same condi- 
tions. They used a quartz-fibre anemometer to measure the velocities 
and manganese-constantan thermocouples to measure the temperatures across 
the boundary layer. 

Pohlhausen [1] appears to have been the first to suggest a 
similarity transformation of the boundary-layer equations. He determined 


the unknown conditions at the surface from the experimental results of 


The numbers in [ ] denote the references. 
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“ Sennide and Beckmann and obtained a series solution for the transformed 
boundary-layer equations. However, his results are restricted to the 
flow conditions applicable to Schmidt and Beckmann's experiments. To 
overcome this restriction Saunders [1] suggested using polynomial approx- 
imations to solve the single fifth-order equation which results from a 
combination of the transformed momentum and energy equations. 

Ostrach [3] presented the first "exact" numerical solutions of 
the transformed free-convection boundary-layer equations for the flow 
about an isothermal vertical plate. He obtained the velocity and temper- 
ature profiles and the heat transfer for a range of Prandtl numbers from 
0.01 to 1000. Ostrach found that his theoretical solutions agreed well 
with the experimental results of Schmidt and Beckmann. Ostrach's work 
was followed by many other numerical solutions encompassing a wide range 
of conditions. A few of the analyses dealing with vertical surfaces are 
briefly discussed below. 

Sparrow and Gregg [4] considered a free-convection boundary- 
layer flow along a vertical plate with a uniform heat flux from the 
surface. They determined a similarity formulation of the equations and 
presented solutions for a range of Prandtl numbers from 0.1 to 100. In 
a later study, Sparrow and Gregg [5] considered non-isothermal vertical 
surfaces which admit to similarity transformations. They presented solu- 
tions for several examples of each of two families of power-law temperature 
distributions. Sparrow and Gregg [6] also presented similarity solutions 
for several cases of variable-property flows. 

Each of the above analyses is concerned with a free-convection 


flow along a surface having a continuous temperature distribution. 
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Shetz and Eichhorn [7] performed experiments using a vertical plate with 
a single step discontinuity in the surface temperature. They used a Mach 
Zehnder interferometer to obtain the temperature and heat transfer data 
for air. They also performed experiments in water using Tellurium dye 
to observe the flow. Hayday, Bowlus and McGraw [8] analysed the problem 
theoretically using the concept of local similarity in eae tipn with a 
transformation similar to that used by Falkner and Skan [9]: this type 
of transformation will henceforth be called a Falkner-Skan transformation. 
Combined free- and forced-convection boundary-layer flows along 
vertical surfaces have also been investigated. Sparrow, Eichhorn and 
Gregg [10] obtained a similarity transformation for this problem and 
presented velocity, temperature and heat transfer results for several 
examples of aiding and opposing flows. Szewczyk [11] approached the problem 
somewhat differently. He set out to determine the effect of free convection 
on forced-convection flows and vice versa. For each case he determined 
the first three terms of a power series expansion where the first term was 
the forced-convection flow in the first case and the free-convection flow 
in the second case. Kubair and Pei [12] again used a similarity formulation 
to determine the solutions for combined free- and forced-convection flows 
of several non-Newtonian fluids. 
Eichhorn [13] analysed the effect of mass transfer on a free- 
convection boundary-layer flow along a vertical, isothermal plate. He 
found that for certain power-law distributions of mass transfer the problem 
could be solved in terms of similarity solutions. Eichhorn considered a 
range of conditions which covered both suction and blowing, and his results 
show that mass transfer has a very strong effect on the heat transfer but 


only a small effect on the skin friction. 
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Most analyses of free-convection boundary-layer flows assume 
that the fluid outside the boundary layer is isothermal. An exception 
is the work of Cheesewright [14] which deals with an isothermal vertical 
surface in a non-isothermal environment. Cheesewright found that there 
are certain temperature distributions of the surrounding fluid for which 
similarity solutions can be obtained. For several examples of these 
temperature distributions, Cheesewright presented the velocity, tempera- 
ture and heat transfer results in graphical form and also tabulated the 
similarity solutions for two examples. 

Each of the above analyses has dealt with free-convection 
boundary-layer flows along vertical plates, but another important surface 
geometry is a vertical cylinder. In analysing vertical-cylinder problems, 
Millsaps and Pohlhausen [15] found that the flow along an isothermal 
vertical cylinder does not admit to a similarity solution. However, 
they did find that for a linear surface-temperature distribution the 
equations admitted to similarity solutions. They presented several numer- 
ical solutions for the velocity and temperature profiles and compared their 
heat transfer results with results obtained by applying integral methods. 
The comparison revealed that integral methods using parabolic approxima- 
tions of the velocity and temperature profiles give heat transfer results 
close to the exact numerical solutions. 

Kuiken [16] also investigated free-convection boundary-layer 
flows along vertical cylinders as part of his analysis of the effect of a 
small radius of curvature on the flow. He considered the flow past thin 
vertical cylinders and slender vertical cones. He aepited a Falkner-Skan 


transformation and used a series expansion to obtain solutions for several 


non-linear surface-temperature distributions. 
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Braun, Ostrach and Heighway [17] extended the analysis of free- 
convection flows along vertical bodies to a wider class of two-dimensional 
and axisymmetric body shapes for which similarity solutions could be 
obtained. Using an integral method, they determined the growth of the 
boundary-layer thickness and the heat transfer for several specific body 
shapes for isothermal surfaces. Non-isothermal surface-temperature distri- 
butions for which Sandi até Jaren netionie could be obtained were also discussed. 

One of the earliest theoretical investigations of free-convection 
boundary-layer flows along a horizontal surface was the analysis of 
Stewartson [18]. He presented a similarity formulation of the equations 
and also obtained a solution for an isothermal surface. However, in his 
analysis of the equations, Stewartson concluded that a boundary-layer flow 
on an upward-facing, heated, semi-infinite horizontal surface was impossible 
but such a flow was possible on a downward-facing, heated, semi-infinite 
horizontal surface. These conclusions were the reverse of what Stewartson 
should have concluded, and this error was pointed out by Gill, Zeh and 
del Casal [19] who showed that the sign of one term in Stewartson's analysis 
was incorrect. 

Rotem and Claassen [20] further examined boundary-layer flows 
along horizontal surfaces by considering similarity formulations for power- 
law surface-temperature distributions and determining the limits of the 
exponents for which solutions canbe obtained. They specifically discussed 
the isothermal and uniform-heat-flux temperature distributions, but they 


presented numerical solutions for the isothermal case only. They also 


considered solutions which are independent of the Prandtl number provided 
that the Prandtl number is tending to the appropriate limit. Im addition 


to their theoretical analysis, Rotem and Claassen performed some experiments 
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using semi-focusing colour-Schlieren photography to obtain temperature 
data and what appeared to be reasonable evidence for the existence of a 
boundary-layer flow on a horizontal semi-infinite surface. They used a 
variety of plate sizes and covered a range of Rayleigh numbers from 0 to 
| 40,000. In a second paper, Rotem and Claassen [21] extended the above 
theoretical analysis to axially symmetric flows. 

The vertical and horizontal surfaces are rather special cases 
of a more general class of inclined surfaces. It appears that the earliest 
investigation of free-convection flows along an inclined surface was that 
of Rich [22]. Rich suggested that the boundary-layer equations for a 
flow along a vertical surface could be modified for inclined surfaces 
simply by modifying rhe buoyancy term. He made no attempt to include the 
effects of pressure variations, and he made no attempt to solve the 
equations. Instead, Rich examined the flows experimentally using a Mach 
Zehnder interferometer to determine the heat transfer from flat plates 
heated isothermally and inclined at angles up to forty degrees from the 
vertical. 

Levy [23] was perhaps the first investigator to present some 
theoretical solutions to the problem of free-convection boundary-layer 
flows along inclined surfaces. He was concerned with the applicability 
of integral methods to free-convection problems, and one of the examples 
he referred to was the inclined isothermal plate. He concluded that the 
accuracy of the heat transfer results was good, and he found his results 
for the inclined plates were in agreement with the experimental results 
of Rich [22]. Michiyoshi [24] also applied an integral method to a free- 


convection flow along an inclined surface. He considered the heat transfer 


from an isothermal, infinitely-wide, thin plate of elliptical cross section. 
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In his analysis, Michiyoshi appears to have neglected the effect of the 
longitudinal pressure gradient and, therefore, his results cannot be 
expected to be accurate for inclinations near the horizontal. 

Kierkus [25] used another approach to the solution of the 
boundary-layer equations for free-convection flows along inclined sur- 
faces, namely a perturbation analysis. Using the similarity solution 
for an isothermal vertical plate as the zeroth-order approximation and 
using a small parameter relating the surface inclination and the Grashof 
number as the perturbation parameter, Kierkus determined the first-order 
approximations for angles up to forty-five degrees on either side of the 
vertical. To check the validity of his analysis, Kierkus also experi- 
mentally determined the velocity, temperature and heat transfer results. 
He found good agreement between theory and experiment and concluded that 
his analysis was valid. 

Another more recent experimental analysis of free-convection 
flows along inclined surfaces was conducted by Hassan and Mohamed [26] 
who were primarily interested in the heat transfer results. To obtain 
these results, Hassan and Mohamed used Boelter-Schmidt heat flux meters. 
Their experiments were performed in air for isothermal surfaces and the 
results cover a range of surface inclinations from 0° to 180°. Wherever 
comparisons were possible, their results were in good agreement with 
other experimental results [2, 22, 25] and with the available theoretical 
résuiltts (3, 25)% 

1.2 Previous Analyses of the Stability of Free-Convection Boundary-Layer 


Flows 
The above discussion indicates that the interest in free- 


convection flows is considerable. The analyses considered above were all 


x 


ods) 20, sasha sid) basookzen me 
od stoma satbeort 8A , 270193 


-lesnestred ofa wasn” enoks | 
ote io. hokteloa Ste: 02) inaoeelel sortons peak! 

“sie bontlank gaols awol? yee: <A -h res ai Salat 
notsuloe vitrsLimte snd goied -sbeylenk un bane 

hae moktemkxorqgs rebr0“t ures oft en sTslq Ips = 
todaexD odd bos wotsentiont soe ide Liat sankrolsy <j dadahone dbeelen 
wbte-se7i72 ons bentarssob eunrsiz: add Smeieg detuned f 
silt Yo obie costs no eagageb Pa ere Ot qu ao ign® x02 anak dis sith 
~tr19qxs' coals ails ly .siaylens aid Re warbi kes sd2 gens ‘oT fe > 
.etivess teténst3 taod be stujsiagne4 eviftsofsy soit Banhintstob® eid ‘ 
tends bobulonoo bas jpemtisq4s bas yrosts Tt ee 4 ” ¥ a 
 biisv eaw's 
mplLIpevsoo~sest to etevisos I[sjasmt L9GXS 31999 eee en Ny “a 
fdS] 5eaedoM bas taesal vd basaubnos ABW sedeiitum neteienii 
nissdo oT .atiueox totedsyd des sda a boaesaasnb ettaatsa Hs 
atevom xuh? dssd $Pkmloe~res Laat beau bpd baw caisaist -sdlseer si 
ada bos asokixve Lamrsiitoet tot rks oF bomltaeg angen - a 
weveusdW) .°O8L of °O moxt ‘enolisniiont soatee 20a 8 evo |e asi 
fitiw Iasuesxgo boog nkversw adivgst viedy yeldkesta axel a do 
isotsetosds sldelinve ala dakw bos sie daa es seh vier srasigpa 


RS sea ch 
Li : 0 


Sry ae ak, rent a og. 


one coe 


concerned with steady, laminar flows; however, another important aspect 
of free-convection flows is the stability of the flow. Some of the 
analyses of stability and transition to turbulent flows are discussed in 
this section. 

Saunders [1, 27] was one of the earliest to study transition, 
but he pointed out that in 1922 Griffiths and Davis had noted the trans- 
ition to turbulent flow on a heated vertical cylinder. Saunders [27] 
performed his experiments on a heated vertical plate placed inside a 
pressurized tank. He was primarily interested in determining the effect 
of pressure on the heat transfer and his study of the transition to a 
turbulent flow was only a minor aspect of his work. However, this study 
did reveal that the position of the transition is very sensitive to 
draughts. Since his apparatus was relatively free from draughts, Saunders 
obtained a higher value for the transition Grashof numbers than that 
obtained by Griffiths and Davis. 

Eckert and Soehnghen [28] were among the earliest to perform 
experiments with the primary purpose of studying the stability and the 
transition to turbulence of free-convection boundary-layer flows. They 
used a Mach Zehnder interferometer to study the flow along an isothermal, 
vertical plate. They observed the formation and amplification of a wave 
disturbance which was initially two-dimensional and purely sinusoidal. 
However, as the disturbance amplified further, the wave form became more 
complex and eventually gave way to a three-dimensional disturbance. On 
the basis of their observations, Eckert and Soehnghen concluded that the 
flow initially becomes unstable due to some small wave disturbance having 


a particular wavenumber. 
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The first theoretical attempt to predict the stability of free- 
convection boundary-layer flows was made by Plapp [29]. Using the eas 
disturbance theory, Plapp derived the disturbance equations including the 
coupling between the disturbance momentum and energy equations. He 
accounted for the effect of surface inclination and demonstrated how to 
account for variable properties. However, since he lacked the necessary 
computer facilities, Plapp could only obtain approximate solutions of the 
uncoupled equations for the case of constant-property flows. In particular, 
he obtained the neutral stability curve for the uncoupled equations for a 
free-convection flow over an isothermal, vertical plate. 

Birch [30] performed another experimental analysis of a free- 
convection flow about an isothermal, vertical plate. He studied the 
stability of the flow by introducing small disturbances at various dis- 
tances from the leading edge of the plate. He introduced these distur- 
bances by passing an electrical pulse through a wire placed in the flow, 
and he observed the flow with a Mach Zehnder interferometer. 

Another experimental technique was utilized by Eckert, Hartnett 
and Irvine [31] for the purpose of detecting three-dimensional effects 
in the transition process. Interferometers cannot detect these effects 
Since the density of the observed fluid is integrated along the light 
beam. Therefore, Eckert, Hartnett and Irvine used smoke threads to 
obtain information to supplement the interferometric studies of free- 
convection flow about an isothermal, vertical plate. They also used this 
technique to determine the effect on the stability of the flow of placing 
an obstacle in the flow at certain positions along the surface. They 
found that the critical Grashof number decreased as the distance of the 


obstacle from the leading edge increased. 
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A theoretical and experimental analysis of the stability of a 
free-convection boundary-layer flow of water around an isothermal, vertical 
plate was presented by Szewczyk [32]. Szewczyk used expansions about the 
critical points to obtain solutions to the uncoupled disturbance equations. 
For his experimental work, he used a thermocouple probe to measure the 
temperature, and he used dye-injection to observe the flow field. 

Nachtsheim [33] was the first to obtain a numerical solution of 
the coupled disturbance equations. Previous investigators had always 
assumed that the effect of coupling was small, but Nachtsheim showed that 
this effect was very significant. He determined both the coupled and 
uncoupled neutral stability curves for free-convection flows of air and 
water along an isothermal, vertical plate. The coupled results revealed 
the existence of two modes of instability: one in which energy is trans- 
ferred to the disturbance by Reynolds stresses, and one in which energy 
is transferred to the disturbance by the interaction of the buoyancy 
forces with the velocity fluctuations. 

The disturbance equations for a vertical plate having a uniform 
heat flux from the surface were derived by Polymeropoulos and Gebhart [34]. 
Polymeropoulos and Gebhart solved the uncoupled equations and showed how 
to convert the results for an isothermal plate for comparison with the 
uniform-heat-flux results. An experimental analysis of this problem was 
also performed by Polymeropoulos and Gebhart [35]. They introduced 
artificial disturbances by means of a vibrating ribbon placed in the flow, 
and they determined the approximate position of the neutral stability 
curve by carefully observing the disturbances with a Mach Zehnder inter- 
ferometer. The results were in good agreement with Nachtsheim's theoret- 


ical results [33] and definitely revealed the importance of coupling. 
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Knowles and Gebhart [36] solved the disturbance equations for 
the case of a free-convection boundary-layer flow along a vertical surface 
having a uniform heat flux. They showed that a thermal-capacity coupling 
exists between the fluid and the surface and that this coupling has a 
first-order effect on the solution. They solved the equations for several 
values of the thermal capacity for air, and the results showed the impor- 
tance of this coupling. The results were also in good agreement with the 
experimental results of Polymeropoulos and Gebhart [35]. In addition, 
Knowles and Gebhart chose a zero-thermal-capacity surface and obtained 
EBiicdtont for sefeieall Prandtl numbers. 

In another analysis of the stability of a free-convection 
boundary-layer flow along a uniform-heat-flux vertical plate, Dring and 
Gebhart [37] considered the spatial amplification characteristics of the 
disturbance as it moved downstream. They determined the constant-amplifi- 
cation-rate contours from which the relative amplification was determined. 
Then using the amplitude of the neutrally stable disturbances as unity, 
Dring and Gebhart presented the amplitude ratio contours. Their results 
revealed that low-frequency, long-wavelength disturbances begin to amplify 
first but higher-frequency, shorter-wavelength disturbances amplify much 
faster. To check these theoretical predictions, Dring and Gebhart [38] 
conducted an experimental analysis. They used a hot-wire anemometer to 
obtain the amplitude and phase profiles of the disturbance velocity and 
to measure the spatial amplification rate as a function of the frequency. 
To measure the amplitude and phase profiles of the disturbance temperature, 
Dring and Gebhart used an interferometer. The relative spatial amplifi- 
cation of the disturbance temperature was determined from interferometric 


moiré patterns. The experimental results were found to be in good agreement 
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with the theoretical results, and they verified that the disturbances 
which are amplified fastest have much higher frequencies than the distur- 
bances which begin to amplify first. 

Hieber and Gebhart [39] also considered the stability of free- 
convection boundary-layer flows along a uniform-heat-flux vertical plate. 
Using a simplified numerical procedure, they were able to extend the 
analysis to a much larger range of Grashof numbers. For a Prandtl number 
of 0.733 as the Grashof number tends to infinity, Hieber and Gebhart found 
that the effect of temperature coupling vanished more rapidly than the 
effect of viscosity. They also noted that the upper branch of the neutral 
stability curve was oscillatory but tending to a non-zero asymptote 
typical of an inviscid instability. For higher Prandtl numbers, the two 
instability modes merged to form a loop in the neutral stability curve. 

As the Prandtl number tended to infinity, the temperature coupling was 
dominant in the instability. For small Prandtl numbers, there was only 

one mode of instability with the temperature coupling effect being negli- 
gible or highly destabilizing for large or small thermal-capacity walls, 
respectively. They also presented empirical correlations between the 

small disturbance theory and the regions in which the flow first becomes 
significantly oscillatory and in which the flow first departs significantly 
from a laminar flow. 

Hieber and Gebhart [40] further examined the above problem in 
an attempt to determine the physical effects associated with the two 
instability modes as the Prandtl number tends to infinity. They used 
asymptotic expansions of both the mean-flow equations and the disturbance- 
flow equations to accomplish their goal. From their analysis, Hieber and 


Gebhart noted that while the boundary-layer theory predicts the formation 
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of an inner or thermal layer and an outer or viscous layer, the stability 
theory combined with the empirical correlation with the experimental 
results predicts that the flow will be turbulent before the two-layer 
structure can fully develop. Hieber and Gebhart also converted the 
results to the isothermal case. 

All of the above analyses have dealt with vertical surfaces, but 
some experimental work has been done on the stability of free-convection 
flows bidet inetinad surfaces. Tritton [41] ees a fibre anemometer to 
study the stability and transition of a free-convection flow along an 
isothermal plate inclined at angles up to fifty degrees on either side 
of the vertical. He noted that the buoyancy effect stabilizes the flow 
for downward-facing surfaces and destabilizes the flow for upward-facing 
surfaces. The experimental results of Lock, Gort and Pond [42] clearly 
showed this effect. A Schlieren apparatus and thermocouples were used 
by Lock, Gort and Pond to obtain a plot of the critical Rayleigh number 
versus the surface inclination. 

Lloyd and Sparrow [43] later confirmed the results of Lock, Gort 
and Pond; but in addition, they revealed a fundamental change in the 
mechanism of the instability as the surface inclination passes through a 
range between 14° and 17° from the vertical for upward-facing surfaces. 
They found that for angles less than 14° the instability was due to wave 
disturbances but for angles greater than 17° the instability took the form 
of longitudinal roll vortices. To study the stability of the flow, Lloyd 
and Sparrow used the same technique used by Sparrow and Husar [44] in a 
study of the formation of roll vortices, namely an electrochemical flow 


visualization technique. 
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1.3 Present Analysis 


The above discussions of steady or mean, laminar, free-convection 
boundary-layer flows indicates that the theoretical solutions for vertical 
(90°) and horizontal (0°) surfaces are well established under a wide range 
of conditions. For other inclinations, solutions are available only by 
approximate integral methods or perturbation analyses under limited condi- 
tions. It is felt that the profiles obtained by integral methods may not 
be sufficiently accurate for the stability analysis, and the perturbation 
analysis has a limited range of application. Therefore, one purpose of 
the present analysis is to attempt to find a boundary-layer formulation 
of the equations which is valid for two-dimensional and axisymmetric flows 
of any Newtonian fluid. along surfaces having smooth, continuous temperature 
distributions. (A finite number of finite discontinuities in the surface 
temperature might be considered by a method used by Hayday, Bowlus and 
McGraw [8].) Before attempting the boundary-layer formulation, it is also 
necessary to establish a normalization procedure and an order-of-magnitude 
analysis which are valid under the conditions stated. A procedure for 
solving these boundary-layer flows must also be found. In particular, 
numerical solutions are sought for constant-property, free-convection 
boundary-layer flows of air along long, inclined, isothermal, plane sur- 
faces. The resulting temperature, pressure, velocity and heat transfer 
data obtained from this analysis will be compared with available theoret- 
ical and experimental results. 

A second purpose of this analysis is to study the stability of 
free-convection boundary-layer flows subjected to two forms of small 
disturbances: a two-dimensional wave disturbance, and a set of longitu- 


dinal roll-vortex disturbances. For the somewhat general mean-flow 
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conditions stated above, an attempt is made to establish the linearized 
forms of the disturbance equations based on the small disturbance theory. 
As in the case of the mean-flow equations, a normalization procedure and 
an order-of-magnitude analysis have to be established such that they are 
valid under the conditions stated for the mean flow. Then procedures for 
solving these equations must be chosen. Solutions are then sought for 
the stability of constant-property, free-convection boundary-layer flows 
of air along long, inclined, isothermal, plane surfaces subjected to both 
forms of disturbances. In particular, the neutral stability curves are 
to be determined, and from these curves the critical Rayleigh numbers 
will be determined. The critical Rayleigh numbers for both disturbance 
forms will be compared to determine the ranges of surface inclinations 
over which each disturbance form governs the stability of the flow. 
Finally, the critical Rayleigh numbers will be compared with the available 


theoretical and experimental data. 
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CHAPTER Il 
FORMULATION OF THE PROBLEM 


DARN Steady-Flow Equations 


Steady, two-dimensional or axisymmetric free-convection flow of a 
Newtonian fluid about. a heated surface or in a buoyant jet is considered. 
(The flow about any cooled surface inclined at 180° to an identical heated 
surface can be described by the same equations and boundary conditions 
provided that the property variations of the fluid are skew-symmetric 
about the ambient temperature.) An attempt is made to reduce the Navier- 
Stokes equations and the energy equation to a set of equations describing 
a free-convection boundary-layer flow along the surface, or in the jet. 
The equations are to be applicable for a temperature distribution which 
varies continuously with position along the surface, or the axis of the 
jet, or which has a finite number of finite discontinuities. 
2.1-1 Governing peer icns 

For the present analysis, thermal radiation effects are assumed 
to be negligible and the above flow is assumed to be described by the 
following vector equations [45]: 
conservation of mass: 

V. (pV) = 0, 
conservation of momentum: 


pv.WV = -VP + V[(i-2u)V.¥] + (Vu).[(VW) + (WW) I] 
3 
Be Ste > (251>1) 
+ u[V*V + V(V.V)] + 9G, 
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conservation of energy: 


V.VP = V.(kVI) + ((i-2u)(V.V)2 


Cp Vv. VE ee 
‘ oe D 5 


+ uV.{L(VV) + (WW) .].V}-nV. [V2V + (9.9) ] 


In addition to these equations, an equation of state and relations for 
the property variations are required. 

Consider equations 2.1-1 applied to a region such as the one 
illustrated in figure 1. It is assumed that the body force field is due 
to gravitational effects alone. It is also assumed that the pressure can 
be expressed as a sum of the hydrostatic pressure, measured for the case 
of no flow, and a departure from this hydrostatic pressure, that is: 

Berl ae Gi 6 
Similarly, the density is expressed as a sum of the density of the undis- 
turbed fluid far from the surface and a departure, that is: 

P = Poo. + SDs I 


The pressure, lee and the density, 0,4, are related as follows: 


0= -1_ 9Po . 
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where K = K(X) is the curvature in the X direction of the surface Y = 0 


for any Z.. 


The set of equations 2.1-1 can now be written as: 
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The above set of equations is not complete until a sufficient 


number of boundary conditions has been established. However, because of 


the semi-infinite nature of the region, the complete set of boundary 


conditions is not evident. This fact combined with the complexity of the 


equations makes it desirable to seek some simplification to the equations 


with a possibility of obtaining a boundary-layer formulation. Such a 


formulation would eliminate the necessity for specifying those boundary 
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conditions which are not evident. 


2.1-2 Boundary-Layer Equations 


Zt 


In an attempt to simplify the equations, the equations are 


normalized and an order-of-magnitude analysis is performed. This 


procedure 


is given in detail in appendix A. The final result of this procedure is 


a set of equations of the boundary-layer type which may be expressed in 


dimensionless form as follows: 
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2.1-3 Boundary Conditions 


Having posed the problem in terms of a boundary layer formulation, 
it is necessary to establish the applicable boundary conditions. To arrive 
at these conditions, it is assumed that the presence of the heated surface 
does not affect the fluid in the region preceding X = O or in the region 
far from the surface. Furthermore, it is assumed that there is no flow 
parallel to the surface at X = 0. Although experiments indicate that these 
conditions are not exact, there is no detailed information known a priori 
about the flow outside the boundary layer to facilitate the establishment 
of more realistic boundary conditions. However, even if it were possible 


to establish more exact boundary conditions, it must be noted that the 


ee ony 


a 


Coe ag A j CSwttys- oA Cates an 20] On 1 * are 7 


[ (I-*wa) 4 


pt 208A fron pa 20 Ke : s | 


*tettalvaro? yeysl ytebaved s Jo amrsa ak eo ldorg a stoi sobil 


svkwis oT .anotykbaos yisboued: istUssbtage ae senighe yyBens 
Ssesitiue bstsed sd ao saitsesiq oft gods bemvaes. at Sri $tb / 


aokges wild) 10 0 SID cst! a9 ait?! a nin, “ = sos 8 


Axohsg s awotsl pibaeoaee) == sia 
— otly sa j 


23 


Garay Renee ae equations 2.1-3 are not valid near X = 0. Thus by applying 
these equations near X = 0, the solutions obtained will be in error, but 
it is assumed that this error will propagate only a short distance inte 
the region where the boundary-layer approximations are valid. 

The other boundary conditions state that there is no flow across 
the solid surface, there is no slip of the fluid in contact with the 
surface, and the fluid at the surface attains the temperature of the 


surface. Therefore, the relevant nondimensional boundary conditions are: 


x = 0, y>0: yY=l1 


x>0, y = 0: 0 = 0(x) 
u=ve=0 
(2.1-4) 
All x, yr; yY=1 
M=1 
é¢= 1 
K=1 


us-Pa = OF=70 


The condition 0 = 0,(x) for x>0 and y = O implies that the 
surface temperature distribution is known. If the heat flux from the 
surface were known, this condition would be replaced by a condition on 
00 at the surface. 
dy 
2.2 Constant-Property Flow 
2.2-1 Governing Equations 


In general, the fluid properties are functions of temperature 
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and pressure. However, the departure from hydrostatic pressure resulting 
from a free-convection flow is usually of a magnitude which has little or 
ae effect on the fluid properties. If temperature variations are small 
throughout the fluid, the property variations Pesittiie from the temper- 
ature variations can frequently be ignored. 

Consider all fluid properties constant with the exception of 
the density. The density is also considered to be constant in all terms 
except the buoyancy terms, in which the Boussinesq approximation is used, 
that is: 

p-==Pell=BeT=T 3}. 


Using the expression for the density as given in section 2.1-1, it follows 


that: 
Pq a -9,,BO 
The set of equations 2.1-3 now becomes: 
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The corresponding boundary conditions are as follows: 


x=0, y>0: u=P,=0=0 
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xo 05: y = Ox) (2.2-2) 


All x, y> ©: -u=P,=6#2:0 


2.2-2 The Falkner-Skan Formulation 

In order to proceed further with the analysis of the above 
equations, consideration must be given to the geometry of the surface. 
The following analysis is concerned primarily with two-dimensional sur- 
faces. Whenever possible, an attempt is made to outline the procedure 
for asixymmetric bodies and for buoyant jets. 

The equation of conservation of mass, or the continuity equation 
as it is frequently called, as it appears in equations 2.2-1 or 2.1-3 can 
be identically satisfied by the introduction of a stream function, i), 


defined as follows: 


oy _ u and dp _ -(1+Aqy)v, (2.2-3) 
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or ow _ yu and op _ -y(1+Aqy)v, (2.2-4) 
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In addition to the stream function, the Falkner-Skan transfor- 


mation [9] is introduced as follows: 


n+ 3 
5-4w- 
V(x,y) = x7 F(n, €), 
4n+2+8w 
5-4w" 
P (sy) ans I(, Gye (2.2-5) 
and 6(x,y) = x" O(n. 6 )5 
where Sc 
n-2+4u" 
nope . 5-4W 
oF 


and C, is a constant used in scaling the boundary layer thickness in the 

n, € plane. The derivation of this transformation is given in appendix B. 
It is also pointed out in appendix B how the transformation can be modified 
for flows for which it is either impossible or inconvenient to introduce 

a stream function, such as in the case of axisymmetric flows. 


Introducing thetransformations given in equations 2.2-3 and 2.2-5 


into equations 2.2-1 leads to the following set of equations: 
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where the subscripts § and n refer to partial differentiation with respect 
to € and n, respectively. 

Before considering the boundary conditions to apply to equations 
2.2-6, it is important to note that the function, a(x), is chosen such 
that OX) = © (Lb) [1+a(x)] x’, which leads to a very simple transformed 


boundary condition. 
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The boundary conditions now become the following: 
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where €° is some small, but non-zero, value of & used in the lateral 
momentum equation to avoid the singularity which arises in an attempt to 
solve for i at. — = 0: 

Equations 2.2-6 are applicable to a reasonably wide class of 
problems. However, the primary purpose of this work is to consider a 
class of problems for which the curvature of the surface is zero, that is 


q = 0. For this particular class for the special cases of a = O and 


a = 7/2 (with a(x) = 0), equations 2.2-6 are replaceable by ordinary 
differential equations since the transformations given by equations 2.2-5 
can be replaced by similarity transformations. However, since equations 
2.2-6 are valid for a(x) = 0 for a@ = 0 and a = 7/2, it must follow that 

all terms involving derivatives with respect to € are identically zero. 
Similarity solutions were first obtained for a = O and a = 7/2 by Stewartson 
[18] and Ostrach [3], respectively. Similarity solutions for a class of 


problems for which the curvature isMOn-Zero were considered by Braun, 


Ostrach and Heighway [17]. 
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2.2-3 The Lefevre Transformation 

Equations 2.2-6 are applicable to fluids having a Prandtl number 
of the order of unity or higher. It would be most desirable to obtain a 
set of equations which would apply for any Prandtl number, including the 
limits of zero and infinity. Lefevre [46] introduced a transformation, 
hereafter referred to as the Lefevre transformation, which was suitable 
for all Prandtl numbers including zero and infinity. Using Lefevre's 
basic idea, the transformation for the present analysis is derived in 


appendix C and is given by the following: 
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The set of equations 2.2-9 with the boundary conditions 2.2-10 
is only applicable to constant-property flows. However, the Falkner-Skan 
transformation and the Lefevre transformation are not restricted to such 
flows. Given some additional relationships to specify the property vari- 
ations, it would be possible to introduce additional dependent variables 
into the derivations of these transformations, but the resulting transfor- 
mation would not be general since they would depend on the specific 
property-variation relationships. No attempt will be made to give any 
examples for variable property fluids since the present work is primarily 


concerned with constant-property fluids. 


2.3 Disturbance Equations 


The combined hydrodynamic and thermal stability of the steady, 
free-convection fluid flows described in the preceding section is consid- 
ered. Using the method of small disturbances, an attempt is made to 
reduce the time-dependent Navier-Stokes equations and the time-dependent 
energy equation to a set of disturbance equations which determines the 


stability of the flow. 


2.3-1 Governing Equations 


The small disturbance theory assumes that each of the dependent 
variables is expressible as a sum of a steady-state component and a small 
disturbance component. In the present analysis, the steady-state (or mean) 
component is indicated by a bar over the variable and the disturbance 
component is indicated by an asterisk. In formulating the analysis, the 
disturbance components are assumed to be three-dimensional and time- 
dependent. It is also assumed that the disturbance components are suffic- 


iently small to permit the neglect of terms which are quadratic or higher- 
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order in the disturbance components. Therefore, the governing equations 
are the three-dimensional forms of equations 2.1-1 with the addition of 
the unsteady-flow terms. Assuming that the mean components satisfy 


equations 2.1-2, the resulting set of equations is the following: 
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The semi-infinite extent of the region in the X direction and 
the lack of distinct boundaries in the Z direction present difficulties 
in attempting to establish a set of boundary conditions for the above set 
of equations. However, even supposing that the boundary conditions are 
available, the complexity of the above equations makes any solution 
attempt very difficult. Therefore, it is desirable to seek some simplifi- 


cation to the above equations. 


2.3-2 Simplified Disturbance Equations 


In appendix D, the set of equations 2.3-1 is normalized and an 
order-of~magnitude analysis is applied. The end result is the set of 
equations D-1l which constitutes one form of the simplified disturbance 
equations. The mean-flow quantities appearing in the equations are assumed 
to be known. However, the equations cannot be solved without having 
auxiliary relationships for M*, c* and K* in terms of the other disturbance 
variables. 

Assuming that these auxiliary relationships are available, then 
a suitable representation must be found for the time dependence, for the 
semi-infinite extent of the region in the X direction and for the lack of 
boundaries in the Z direction. This problem is approached by assuming 
that the disturbance variables are proportional to exp [i(Ax* + Qz - BT)]. 


Using this assumed proportionality, the disturbance variables are bounded 
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as X or Z goes to infinity, and the stability of the flow can be analysed 
by considering a single component of a Fourier series in x* and z. The 
real parts of Xi and Q, d- and Qs are the wave numbers in the x* and z 
directions, respectively, and the real part of B, Bis is the frequency of 
the disturbance. The imaginary parts of A, 2 and B, hie 2. and Bi» are 
amplification factors. The disturbance grows with increasing x*, z and T 
if, and a. are negative and B. is positive. If ye Q. and B. are all 
zero, the disturbance neither amplifies or attenuates and the flow is said 
to be neutrally stable. 

For the stability analyses of free-convection flows, Knowles and 
Gebhart [36] have shown that for a = 1/2 free-convection flows subjected 
to two-dimensional disturbances have a lower critical Rayleigh number than 
the same flows subjected to three-dimensional disturbances. However, 
since the present analysis is for all a, an attempt is made to solve the 
disturbance equations as they apply for two types of disturbances. The 
first type is a two-dimensional disturbance in the form of a wave which 
travels in the X direction with an amplitude dependent only on y. The 
second type is a three-dimensional disturbance in the form of a set of 
roll vortices having their axes parallel to the X direction. These vortices 
are periodic in z and have amplitudes which depend only on y. In addition, 
these vortices are assumed to be independent of x*, but this does not imply 


that u* is zero. 


2.3-3 Tollmien-Schlichting Wave Disturbances 


In forced-convection stability analyses, two-dimensional wave 
disturbances similar to those described above are called Tollmien-Schlichting 


waves. This terminology is also applied in the present stability analysis. 
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For Tollmien-Schlichting wave disturbances, the set of equations D-1l, with 


the chosen function for H(W), can be simplified to the following: 
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(ity ,)c(1+a) [57 + (itAqy) ox* + by! rg { 'd UGeagy). ag cul 


ah) 4 V(1ta)90* | uk(1+a) 96, | y,*c, u(lta) 30 ud da 
jae bes oy a (1+Aqy) 5x eo CisAqy) Ox (1+Aqy) dx 


+ 


v(1+a) 96, i (1+74)c*, u(1+a) 90 ud da v(1+a)36)4 
oy (1+Aqy) 0x (l1+Aqy) dx oy 


1 2 
-5 (1+w~) 25 5 5 SNe 
(lta) Ra. ° * -K* 0°80 Aq 30 Jsina 30. K,0°0% 
L {lage Cuekay). os) Cidiyaiaa) oy) 892 
fence ae, 208 i a79* .  Jsina _90*, | dK*98 | OK 96% 
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Assuming that the disturbance variables are proportional to 
exp[i(Ax*-BT)] andthat the stability of the flow can be analysed by consid- 
ering a single component of a Fourier series in x*, the disturbance variables 


are expressed as follows: 


i(Ax* - BT) 


u*(x*, y, T) 6, (y) e 


vk (x%, Y> T) = 6, (y) see’ ae 


i(Ax* - BT) 


O*(x*, ys TF s(y) e 


eat yer) @ OF eP 


MRC, y, T) = AG) ero 7 OY 


i(Ax* - BT) 


c*(x*, y, T) = Aty)re 


i(Ax* - BT) 


K¥(x*, y, T) = Ty) e 


Substituting these expressions into the set of equations 2.3-2 and 
eliminating the exponential factor from each term, the following set of 


equations is obtained: 
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1 2 
See: eS 2 
Ray shee u Guy Jycosada, , es 
(jtJysina) ‘*(1+Aqy) ‘dx ~ dk’ = 
-§ (1+?) = ide Sy Glesie) = 
a Ra [ uf ce - AQV, 4 ae + Ra, : var 0 
(1+Aqy) ‘dx Oy dy 
ma rae me 2 % he weds ra 
oRa, 5 (Ho ia 36 Ra pt irw") 200M 20MAq (1+y,)v oMJsina ,d“6 1 
Pe =< d raga mae tis or pay 
dy Oy (1+Aqy) (jtJysina)* dy 
i(l+y,),B Au Ra ~$(140") 037M . 200M OAq OM 
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So) ee, i es — 
~ (jtJysina) dy (1+Aqy) 7 (1+Aqy)“ ~~ (jtJysina) (1+Aqy) 


_ MJ" sin®*a _ (1+7,)av _ (1+¥q) (2Aqv Fe du : v8¥q, }46] 
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{ (1+¥q) ;iAqB _ iddu _ 2iAAqu rae idVa)B omen 
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Ra 75 (1+u*) oAq 32M o 202 + A2q2. aM 
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Se — Sa = 
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Giciasy, oe Cag) Oyo as as 
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ah A eee 
: 


46 


if. ore. z 
_ 2(1+Yg)A qyudg , __1 ata ou Aqvy 1} 9) 


(1+Aqy)” dx (l1+Aqy) dy 


peel 2 w pa ee 
ioR 5 (1+w") peg = 2 g Ra, 5 (1+w*) __ioMAdg ioMAJsina 
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2 i (47g) Av, — (+Y¥g)Aqu _ (1474) 3u; 462 4 ¢+yg) 1. TeX 
(1+Aqy) * (1+Aqy) dy dy (itAqy) ‘ 
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1 2 
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1 2 

= = 2h(14w?) .- = 
2qdu ie A*qyu? dq,; r , ORa 5 du Aqu 

+ Cisse) diag)? age tt caame oe 


ou | ya 
dy  (1+Aqy) dy” 


ae 1+) - es ma 
oRa 5 ( dA.297u A du A7qu 
+ < + 1. SE 
ay! ay? * (tAqy) dy ~ (i+Aqy)7! 
a 2 f. a 2 fs 
, ORa 5 (1+w") Jsina dA du Aqu ye oRa 5 (1+w a 2° 
2 (j+Jysina) dy dy ° (1+Aqy) oy? 
07u A*q? du A3gq3u oRa -3 (14+w*) AJsina ,07u 
(itAqy) dy? ~ GtAqy)? dy * G+aqy)?)* + Gitysina) ‘ay? 
2Aq uy 2 ORa, 5 (1+w" ) AJ*sin2a peo i" Aqu } 
(1+Aqy) dy (jt+Jysina) ” oy (1+Aqy) 
-§(1+w*) 42 = = 
4 ORa . ASA pours Aquiey = 0 
(1+Aqy) “dy  (1+Aqy) 
(2.3-3) 
ie Cay ee ; 
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dy” (1+Aqy) ~ (jtJysina) * dy 
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1 2 atiez i d 
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ud_da , V(1+a)90) , Ray re (14a) T 3% esta 00) 
(1+Aqy) dx oy Oy” " (L+Aqy) dy 


Jsina 390 Ra, 5 (lw )(14a)96 aT _ 


Spottt ee 
~ (jtJysina) By! 45 dy dy 0 


The above set of equations can be employed in the stability 
analysis of any free-convection flow which satisfies the boundary-layer 
equations 2.1-3 and the boundary conditions 2.1-4, assuming that the distur- 
bance has the form of Tollmien-Schlichting waves. However, auxiliary 
relationships for r, A, A and T and a set of boundary conditions must be 
established before a solution can be attempted. 

By considering constant-property flows over two-dimensional 
surfaces, equations 2.3-3 can be replaced by a set of more simplified 


equations. From the Boussinesq approximation, 


p = 0,.(1-B9). 
However, since p =p + px 
and © = 6 + O*, 
then P = p,,(1-B0) 
and p*. = 7-4 pi BO*. 


Using the normalized forms of these expressions along with the above 
assumptions and assuming e<<l where € = BO. the disturbance continuity 
equation reduces to a form which can be satisfied by introducing a distur- 
bance stream function, ~*, defined by: 
ux = dp 
oy 
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Utilizing the above assumptions and introducing the disturbance stream 


function, equations 2.3-2 reduce to the following: 


5 (1+? da 


tl % AyRay,__ dq daye! 0? pe s Aq. 0*W 
Sees 297 (1+Aqy) dx er = OT acta dyoT 


(1+w? 2 
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1 2 
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————— ee + 
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1 (14w? f 3h) 
= ORA-* A3q? A ae 2Aq__ 9 
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3A7q* a2 2Aq 9 Spx i atye 2 9 tye 
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pte aa 2 
fs any) + oRa, : jAgsina 6* | Ra) Apa sina dav* 
dy (1+Aqy) ) (1+Aqy) dx 
in 1 09% cosa , 00% sino, 
(1+Aqy) ox* oy 
(2.3-4) 
- : 5 (Luar 
(1+a)00* (1l+a)u 06% e “Cray aod ce 7 Ra ud* da 
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—$(1+w?) 2 2 
a Ra, (1+a) 1 0“O% a“O% Aq d0%* 


Z (itAqy)? oxe™ * dy2 * CitAqy) dy | 


Following the procedure used to obtain equations 2.3-1, it is assumed that 
W* and 8* are proportional to exp [i(Ax* - Bt)] and a single component of 
a Fourier series in x* is used to analyse the stability of the flow. 


Therefore, W* and 6* are assumed to have the following forms: 


ot (Ax* - BT) 


We (x*, y, T) oy) 


i(Ax* - BT) 


and Oke, ya) 0) = Sy) ae 


Substituting these expressions into equations 2.3-4 and eliminating the 


exponential factor in each term, the set of disturbance equations becomes: 
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2 we 
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al 2 Be. si 2 cs 
(lta) B _ Au es iRa, 5 (rue) u da) s s + iRa, 5 (iw ) (ita) vas 
(1+Aqy) (1+Aqy) a dy 
a —§ (1+w7) 
A(i+a) 96 > , iRa (1+a) 90 + Oda,do 
(1+Aqy) dy (1+Aqy) ox dx° dy 
oe | 2 
_ ina, 8") 14a) aqeee ce ~ Fyyterey 6} 
dy” " (1+Aqy) dy (1+Aqy)? 


It is possible to simplify the disturbance equations further by 


assuming that the mean flow is parallel. 


This assumption implies that the 


mean-flow velocities and temperature do not vary with x and that v is 


approximately zero. 


Utilizing this assumption, equations 2.3-5 reduce to: 
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Aq Aq? + nya — 3A2q? 
P (1+Aqy) ( recone ) 9} 
(2.3-6) 
eel | 2 x Fares 2 

j (ta) 8 - AU ie , iuRay 5 (1-+w Yaa ar A(1t+a) 38> iORay, 5 (1+ Jaa dd 
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Before attempting to solve either equations 2.3-5 or 2.3-6, a set 
of boundary conditions must be established. Firstly, from the assumption 
of no flow along or through the surface the disturbance velocities vanish 
along the surface y = O and far from the surface. Secondly, the distur- 
bance temperature is assumed to vanish far from the surface. These 


conditions can be stated as follows: 


(223-7) 
and yro 2 @ =do=s =0 


The remaining condition on s is not easily specified. As Knowles and 
Gebhart [36] pointed out, the second condition for s depends on the thermal 
capacity and the thermal conductance of the material used for the heated 
surface. If the thermal capacity of the surface material and the thermal 
conductance normal to the surface are large, it is possible to maintain an 
isothermal surface, implying that the disturbance temperature is zero at 
the surface. If the thermal capacity of the surface material is zero and 


the thermal conductance is negliglible in the direction parallel to y = 0, 
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it is possible to maintain a uniform heat flux along the surface, implying 
that the disturbance heat flux is zero at the surface. The inpermed(are 
conditions can be obtained by considering an energy balance on a small 
element of the material between the heaters, used for heating the surface, 
and the boundary, y = 0. Consider the energy balance after subtracting 
the mean-flow energy terms and neglecting the disturbance conduction terms 


in the direction parallel to y = 0, that is, 


b C 30% k 30% 
pe | TS 


0 
Otte oY =0 
2 2 
or s(o) = iL Ra, aes )ds(o), ee 
Bb dy (2.3-8) 


where b is the thickness of the surface material and it is assumed to be 
small relative to L. The above conditions 2.3-7 and 2.3-8 are applicable 


to either equations 2.3-5 or 2.3-6. 


2.3-4 Taylor-Goertler Roll-Vortex Disturbances 


The roll-vortex disturbances described in section 2.3-2 are 
called Taylor-Goertler roll vortices as are their counterparts in forced- 
convection flows. These roll vortices are independent of x*, and there- 


fore, the set of equations D-1 can be simplified to the following: 


oats (1+Y 4) ,ov* Aque 1 Ow 4 v*Jsind) | 4 v*dVq 
OT dy (1+Aqy) (jtJysina) ‘dz dy 
RIO) ey ae 
4 Ra, +¥q*vJsino ui (ai mn Jycosada, -(1+¥q)us 


(jtJysina) (jtJysina) ‘dx ~ dx’ © (1+Aqy) 


Teme B m0 sompled Yavene ns sntsebtenee xd oe od io 
.sotiwe siz gntssed rot been .exo spor silt apauged a aia oi 3 
 gabsoesidye duet Suaiial wgyeds ond yebienod 0s et “ ibenal ds 
emred notioubhoo soasdiwielh oft gniiaelgen bap ectiss yarens worsen 

eh aed .0 = y oF Lolleseg nols soak‘ 3 


ol 8 : vel 8 od dq 


OQ = 36 47 


= 


| : | hie : 
ee ae) eee ed ee 
ojeb pa id: = (ole 
(8+€.9) eb 4 add 


4 


a 


od oJ bemges ab 32 bos Iebrstsi.ssaivwe 919 Yo Besndoldy od, a dom 
sidepiiggs ois 6-€.8 ban T42.8 ano Hi bties oyods edt 4d og em 
a 

A-E.9 t0 OnE. 8 acolaanps, edits 
i! on 


asonadsusald xst3x0V—-LLoOr wel 


9u6 S~E.S: nokdoee ub badtusesh asonadttath xorwy-Ker bahia 
~beo7e} at ataequogsiuod sheds STA es esolrroy on sal sxepd—todiyat 
-stedd bas .** t0 jopbasqsbat ots e95t4z0v Sto eaee vewol? 4 
tantwollot sia od belRkiqmie 4d na 1-0 anolisups Jo. eins 


55 


Ya%y 1 * du. Aaqv *xOY *Ov  voYa* 
pes St Tan ot qv u%dVq_. Nae CV emmanuel a, 
Citiqyy #88 tay): SyG3," 0 G2? Bax ) > Gigeetnngs ? 7 9 


(1474) ,9*w fee 1 [ 2Aqu du* + Jsinadw _ Uy oYg ow 
dydT (jtJysina)*(1+Aqy) dz ava O20 dy OT 
=2 =s(ltos) ave 5. gaz : 

cs Aqu OYg* bs Ra, “¢t74) (2 w 4 ov. ow 

(1+Aqy) (jtJysina) 92 dy- | dy dy 


# 2 Gt¥g) pt _4d + Jycosada, (wou , udw  _ Aqwu ‘e vo ?v* 
(j+tJysina) * (1+Aqy) ‘dx ~ dk’. oy dy (1+Aqy) ozoy 
+ 2vJsinadw , wuJcosa da. _ Ov dove f wJsinadv, Ss vo¥a Bw) 
dy (l1+Aqy) dx dy oz hay dy dy 
=*(14+07) ,2- 7 one 
Z oRa, . o°M ow 4 1 ous - wisind,, + om, 28 Ww 
Oy dy © (jtJysina) ‘dz dy" dy” 
asain ow 7 _Jsina _ ow a! cAqav* ‘ AqwJsina, 
(1+Aqy) dy (jtJysina) dy (1+Aqy) (j+Jysina) oz \+* 
i 297w , 3Jsinadv* wJ7sin’o, M,0°w Aq 37 w 
+ Geiysinay”‘ 927 * 32 1 + UyyF + Ctaay) ay? 


+ 


A?q?_ dw M _, d°w 2Jsinad?w _ _d°vs Aq 
° (1+Aqyy By! $ (jiJysind) d223y z dy? dy*dz - (i+Aqy) ° 


Q 


Jsinadw 9d*v* A*q? dv* — wJsina M 
at eo al ae es Pi aa ee SE 
‘ oy dydz> 3 (1+Aqy) Oe a long (jtJysina) ( 


co 
sop pea 
a ee 


~ —— = 


rr) = Suny CREE 
w6 mikuyry ae 
ce Ss EM) de I 


26S onketw ~ *y6 t 
aig cK + Se Conemeaaey © 


56 


Jsinad?v*  J*sin*adw . Jsinad2w  d3vx MJ*sin7a wJsina 
2 1Sinoe ses Tsim ye Sy") + ee 
dyoz oy Oz Oz (jtJysina) 


=(07 


_ ove Ocosa Ra oYg* 
OZ as ae E 


85. 


G+¥q) _92ut , du ave | Agu vA) | Ray et ae Ges) (vats 
(jtJysina)°dzdt dy dz (1+Aqy) dz (jtJysina)* dzdy 


1 (au SUF Aqvau* _ 1 __ i , Jycosada, aw, 
(1+Aqy) ‘dx dz Oz (1+Aqy) ‘dx ~ dx’ oT 


1 9¥q ow 1 OYg* -vdu Uae ou F Aqvy)} 
ox 


~ (i+Aqy) ot» (jtJysina) dz By | (tAqy) “ox 


-)(1+w?) 2 = 2 
_ ORa, { 1 3 Me) Bu i Aqu 18% | fe u* 
r (jtJysina) dzdy dy (1+Aqy) (jtJysina) *dzdy 


3 2 22 
¥ Aq du*,0M M d- ux Aq q“u* 1 A*q*du* 
(1+Aqy) 92 ly is (jzaysind) \dzdy" # (1+Aqy) dzdy (1+Aqy ) Oz 


AwJsinadq, , M [+ Jsinad*u* , AqJsina our, 


aig (jtJysina) ” ozdy  (1+Aqy) dz 


I+ 


M d°u* 1 oM* .97u Aq du _——A*q*u 
+ Gitysina)® 2° * Gidysina) dz (ay? 7 Citagy) dy ~ Ctday) = 


2 
- = -w 
Jsina oM* .du Aqu ae Osina Ra, oYg* 


(GtJysina)” ag lsy © (1+Aqy) BO. Oz 


(2,5-9} 


eT 


-(1+y,) Aq dux* 2v*du udv* O7u*  v*d2u 9 dv* du oYg dux* 
q iggy so dy dy dydt Oy H Oy dy oy’ oy (Br 


oe 2 = 
a v*ou i Aquv* 1. nae te ) 1 oY erapout _ 2Aquu* j 
(1+Aqy) (1+Aqy) Eval oT (1+Aqy) 
- (1+Y 4) vey ux a ov du* rf L (2Aqvaus u* d*u i du* du 
dy oy (1+Aqy) dy dyox dy dx 
Aqu* & 2Au* udq | i uu* dq Nn * vo7u 
+ - 
Bp + Casnagye Cas + Gm) ~ Gaeeyy? ag) > TUG 
avout fl gudta 4 du dus Aquov. payee Au 
si dy dy PC Ray) « Oyox oy Ox dy % Ry a atts (Ss i 
2 q3u, A*qyu? - ta vou* u* du , Aqv 
z Poy) ~ (1+Aqy) [ Oy + “(itAqy) Ox a )] 
oY gt vou u , Aqv 
ES oy Lao oy (1+Aqy) GS + I} 


2 = 
_ ona, °C") 92, _Agué ut) _ zak UY G8 Age due 
bs dy? (1tAqy) dy dy- — (itAqy) dy 
A? q?u* 1 OM Aqu*Jsina rn Jsinadu* 
im (1+Aqy) i oe (jtJysina) ay 3 (1+Aqy) — oy 


1 07u* o2Mx du a 207 AB Aq. ou 
(jzdysina) d27 ) * dy? rey dy! ~ dy | by? (1+Aqy) dy 


+ myGy pA 
é } 


< 


(epee aa 


ae 
v6 ca fe 5 


58 


3 A*qu - Jsina 9oM*.du Aqu M,.0%u% 2Aq 07u* 
Civaay)7! (\4@jaayeina) ay loy-aneciragy). cit iloyhal™ Gueaqy) dy? 


a, BS GBS . Sux A*q’u* 15 MJsina ,d*u* | __2Aq dur, 


(1+Aqy)“ dy 5: (1+Aqy) GtJysina) -ay? (1+Aqy) dy 


M J*sin*0du* 03ux* Aq u*J7sin7a, 97ux* 
$ (j+Jysina) oy dyoz- (i+Aqy) © ~ O27 )] 
MJsina 07u*x* M* 92u 2A 07u Raat du A3q?u 
+ < oe at ee ee ee ee 
GtJysina)® O27 [ay + Ctaqy) By? ~ CitAqy)= dy * (tAqy) 


+1 


M*Jsino, (oe 4 2Aq du, M*J?2sin2a (ou 4 Aqu } 
(j#Jysina) dy” © (1+Aqy) dy ~ (jtJysina)7*dy * (1+Aqy) 


—y2 Ad: 2 
4 ORay é ee sino , Yq*Aqsina | Ra, Joa )va*sino do. 
Bo, yey (1+Aqy) (1t+Aqy) dx 


ae an = -1 (14w? he 
(1+¥,)€ (14a) ae ; ne, ona 5 (474) 


_——t da(ud* + u*6) 
(1+Aqy) dx 


u(lta) 96 ud da 


v(1+a)d0* , uk(1ta) 00, | pu(lta) 98 | 
(1+Aqy) 0x (l1+Aqy) dx 


dy * Gray) ox! * | 


+ 


v(1+a) 96, rva*e i (147 4)c*) 


+ 3y 
“E(140) ee 925 : 5 
i (1+a)Ra, “4K (2 x Aq 08 , __Jsina 08 


ay? * GtAay) By * GtJysina) dy! 


59 


K,9*0% Aq -90* , Jsina 1 076% 
2 LSyz SSUthay). oy 4. i LivsinOnoys (jatystnn) 2322. 


4 OK* 98 3K 96%, 
dy dy dy dy 


An approach similar to that employed in section 2.3-3 is used 
in the analysis of the above equations) If it is assumed that the distur- 
bance variables are proportional to exp [i(@z - BT)] and that the stability 
of the flow can be analySed by considering a single component of a Fourier 
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Substituting these expressions into equations 2.3-9 and eliminating the 


exponential factor in each term leads to the following set of equations: 
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Equations 2.3-10 can be utilized in the stability analysis of 
any free-convection flow which satisfies the boundary-layer equations 2.1-3 
and the boundary conditions 2.1-4, assuming that the disturbance takes the 
form of a set of roll vortices which satisfies the assumptions involved in 
obtaining equations 2.3-10. However, a set of auxiliary relationships for 
r, A, A and T and a set of boundary conditions are required before attempt- 
ing a solution. 

As stated in section 2.3-3, the present study is primarily 
concerned with a study of constant-property, free-convection flows over 
two-dimensional surfaces. Using the Boussinesq approximation for the 
density as in section 2.3-3 and assuming that e<<l, the disturbance 
equations reduce to 
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The above set of equations appears to be over-determined since 
there are five equations in four unknowns. However, a closer inspection 
reveals that the fourth equation can be obtained from the second and 
third equations by the procedure which follows. Firstly, multiply the 
third equation by (1+Aqy) and then take the derivative with respect to 
y of the resulting equation. Secondly, take the derivative with respect 
to x* of the second equation, and note that some of tee are of the 
order of the terms previously neglected and, therefore, these terms must 
also be neglected. Finally, add the two modified equations and multiply 
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(1ta)Ra ~801+w7) 42 peas 
- L [Ss + — 4 S - 95] 
dy (1+Aqy) dy 


Following the procedure of section 2.3-3, the mean flow is 
assumed to be parallel. This assumption reduces equations 2.3-12 to the 


following set of equations: 
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1 2 
(1+a)Ra, 5 (14 ) a’ st Aqui daa) Vs, 
dy” " (1+Aqy) dy 


The appropriate boundary conditions to be satisfied by equations 
2.3-13 follow from the conditions that the disturbance velocities vanish 
along the surface y = O and far from the surface. In addition, the distur- 
bance temperature must vanish far from the surface and satisfy the energy 
balance at the surface as given by the expression 2.3-8. These boundary 


conditions may be written as: 


1 2 ash 
2 2 
-5(1+w*) 
and s = iL Ra. ° ds *. 
Bb ay (2.3-14) 
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CHAPTER III 


SOLUTION OF THE PROBLEM 


In chapter 2, the formulation of the problem is presented in 
two steps: the formulation of the steady-flow or mean-flow equations and 
the formulation of the disturbance equations. The resulting sets of 
equations from these steps are distinct in character and, therefore, they 
require different methods of solution. This chapter presents the methods 


of solution adopted for the present analysis. 


3.1 Solution of the Mean-Flow Equations 


There are several possible approaches to solving the mean-flow 
equations. For example, given the property variations, an explicit finite- 
difference scheme might be applied to equations 2.1-3 subject to the 
boundary conditions 2.1-4. An approach such as this was attempted in the 
present investigation but it was abandoned due to difficulties in starting 
the solution at x = O because of the singularity at x = 0 and because of 
the problem of obtaining accurate approximations to the derivatives with 
respect to x in this region. In addition, explicit finite-difference 
schemes have an inherent instability. 

A second possible approach would be to transform the set of 
equations 2.1-3 by a transformation such as the von Mises transformation. 
In the transformed plane, the boundary-layer thickness is more uniform 
throughout the region of interest. In addition to the transformation, 
the explicit finite-difference scheme may be replaced by an implicit 
finite-difference scheme, thereby avoiding any inherent instability. 


However, the von Mises transformation introduces a singularity at the 
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surface, and this singularity may be difficult to handle in the p OMALE Wot 
This approach was also tried, but it was abandoned because of difficulties 
attributed to the singularity. 

A third finite-difference scheme, which is similar to the second 
scheme, is given in a paper by Patankar and Spalding [47]. Patankar and 
Spalding suggested a transformation of the equations of the von Mises type, 
but which uses a new definition of the nondimensional stream function. 

The authors stated that the advantage of their method was that it provided 
a definite bound on the boundary-layer thickness. However, their method 
still has the disadvantage of introducing a singularity at the surface. 

Integral methods offer a fourth possible approach to the 
solution of the mean-flow equations. Levy [23] investigated the possi- 
bility of using integral methods in free-convection problems. Michiyoshi 
[24] also used an integral method in his analysis. The disadvantage of 
integral methods is their inability to accurately predict the shape and 
magnitude of the velocity and temperature profiles. Since the influence 
of these profiles on the stability analysis is not yet clear, it is felt 
that an attempt should be made to keep the errors in the profiles to a 
minimum. 

Another possible solution approach is to use a perturbation 
expansion. Kierkus [25] used a perturbation expansion in terms of a 
small Grashof-number parameter to study the effect of inclination from 
the vertical on laminar, free-convection flow about an isothermal plate. 
Although the results of this perturbation analysis are in good agreement 


with experimental results, this approach is limited by the size of the 


perturbation parameter. 
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The difficulties which may be encountered in utilizing any of 
the above methods suggest that perhaps some other approach might be 
sought. For the transformed boundary-layer equations 2.2-9 subject to the 
boundary-conditions 2.2-10, it is possible to use another approach; however, 
‘the method is not restricted to the flows described by these equations, 
but instead it can be applied to any set of equations expressed in terms 
of a Falkner-Skan transformation. It is noted in chapter 2 and appendix 
B how the Falkner-Skan transformation might be extended to more general 
flows. 
Basically, the scheme chosen replaces the schemes mentioned 
above by a numerical integration, for which there are two broad choices: 
the equations may be integrated as a single set of partial differential 
equations, or the dependent variables can be expanded using, for example, 
a Goertler expansion which leads to an infinite set of ordinary differ- 
ential equations. An important consideration in applying the latter method 
is the number of terms needed in the expansion to obtain the desired 
accuracy in the solution. Although it is believed that either of the 
above methods is suitable, the former procedure is chosen for the present 
analysis. The latter method was used by Kuiken [16] in a study of axi- 
symmetric free-convection boundary-layer flows past vertical cylinders 
and cones. Saville and Churchill [48] also used a Goertler expansion in 
their analysis of horizontal cylinders and vertical axisymmetric bodies. 
The numerical integration in the former method is made possible 
by replacing the partial derivatives with respect to € by finite-difference 
approximations as discussed by Hartree and Womersley [49]. The method 
has been applied to laminar, forced-convection boundary-layer flows for 
some time, but the method has only recently been adopted to free-convection 


boundary-layer problems. After applying a Falkner-Skan transformation, 
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Hayday, Bowlus and McGraw [18] applied this method of solution to'a free- 
convection flow past a vertical plate with a step discontinuity in the 
surface temperature. The apparent success of the method in handling such 
a difficult boundary condition suggested that this approach might also be 
applied to the present work. 

After replacing the partial derivatives with respect to & by 
finite-difference approximations, the equations reduce to ordinary differ- 
ential equations for any value of §. In order to solve the equations for 
a particular €, it is assumed that the equations have been solved at all 
previous steps in the €-direction. At § = 0, it is assumed that the deriva- 
tives with respect to € remain finite, and therefore, these derivatives 
are eliminated because of the coefficient €. However, § cannot be set to 
zero uniformly throughout the equations. This is a consequence of the 
coefficient of the pressure term in the lateral momentum equation. In order 
to start the procedure, € must be given a small non-zero value in this term. 

For the finite-difference approximations to the €-derivatives, a 
two-point difference formula must be used for the first step beyond € = 0. 
For subsequent steps, a three-point formula may be used. The two- and three- 
point formulae were given by Smith and Clutter [50] as well as by Hayday, 
Bowlus and McGraw [8], but they are repeated here for completeness. For 


any variable, say S, the two- and three-point formulae at € = g, are: 
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respectively, where S represents f, f T, X or a. It should be noted 


c? 
that a four-point formula may be used beyond the third step. Smith and 
Clutter experimented with two-, three- and four-point formulae and found 

that the three- and four-point formulae were stable and accurate. In fact, 
they found that the three-point formula was sufficiently accurate such 

that the four-point formula was not required. Since Smith and Clutter 

were concerned with forced-convection flows, it does not necessarily follow 
that the results of these numerical experiments are valid for free-convection 
flows. However, only the two- and three-point formulae are used in the 
present work. 

After suitably discretizing the €-derivatives in the equations, 
the resulting equations may be solved at each € by any of the methods used 
in obtaining similarity solutions. For the present work, the equations 
are forward-integrated, starting at € = 0, using a fourth-order Runge- 
Kutta procedure. To start the integration, it is necessary to guess 
values for fr¢ (0, %&)¢ Xz (0, €) and 7 (0, &), which are the eigenvalues. 
An iterative process similar to that followed by Nachtsheim and Swigert 
[51] is used to determine the eigenvalues such that the asymptotic boundary 
conditions for t > © are satisfied within a specified accuracy, e.g. LO 
After making the initial guesses for the eigenvalues, successive approxi- 
mations are obtained by means of a modified Newton-Raphson procedure which 
incorporates the asymptotic nature of the boundary layer. That is, to 
the conditions 

fr (02 6) = X (os 6) = 7 (Sn, &) = 0, 
the following conditions are added: 

Ere (eos E) = Xz Gor 5) = Ty Geos 6) 2 0; 


where T,, denotes the boundary-layer "edge". A detailed discussion of this 
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procedure was given by Nachtsheim and Swigert, but a brief description 


is also given here. 


and 


After introducing the following definitions 
Pie a frr (0, Sey 


= Xt (0, eH 


a 
NO 
+a) 

| 


By thena(Ogteays 


where €, is the nth step in the €-direction, the first step of the iterative 


process consists of determining 
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Next, this quantity is minimized with respect to the increments AY) 


AT9, and AT3n. The incremental changes in I'j,, Ign and I'3, necessary to 


effect this minimization are then found and used to modify the eigenvalues. 


It must be noted that the above process will converge only if 


the initial approximations to the eigenvalues lie within certain limits 


of the correct eigenvalues, these limits being the boundaries of the 


region of convergence. To facilitate the solution for the eigenvalues, 


the error indicated by Enj is minimized at several 54 <S0o in order to keep 


the set of approximate eigenvalues within a region of convergence. Since 


the value of f, is unknown a priori, this procedure may also be used to 
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find a suitable measure of [. The procedure assumes that for any reason- 
able set of approximate eigenvalues there will be some C1 < So for which 
En] lies within a fairly large region of convergence. It seems reasonable 
to expect that as Cj becomes larger, the region of convergence becomes 
progressively smaller. Having minimized En] and modified the approximate 
eigenvalues, the equations are forward-integrated from zero to i Me 
This procedure is repeated m times until E,m is less than a specified 
value, e.g. HORS. Lhe corresponding value of om then defines €,. Clearly 
this method permits a much wider range of initial approximations to the 
eigenvalues for which the process will converge; however, the choice of 

¢, and the optimum number of steps m can be determined only by trial and 
error methods. 

The above procedure for finding To was used initially; but after 
finding a suitable value, f. was specified to provide a constant value for all 
€ such that for any particular € the necessary information at previous 
values of € would be available for all [ < Gq». The value specified varied 
with the surface inclination and corresponded to a value for yo at §€ = 1.0. 
At a = 0, yo at & = 1.0 was given a value of 10.0, while at a = 7/2, it 
was given a value of 8.0. For other angles, yo at € = 1.0 varied between 
these limits. 

The utilization of the Falkner-Skan formulation in conjunction 
with the local similarity concept has several advantages over other possible 
methods. It provides a definite procedure for starting the solution at 
& = 0, and the boundary-layer thickness remains finite at § = 0 since the 
procedure removes some of the €-variation of the transformed boundary- 
layer thickness. In addition to these advantages, the most important 


advantage is that, for a given €, the equations reduce to ordinary 
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differential equations. 


3.2 Solution of the Disturbance Equations 
3.2-1 Tollmien-Schlichting Wave Disturbances 


This section presents a discussion of some of the numerical 
procedures used in previous work as well as the procedure used in the 
present work. This discussion refers specifically to equations 2.3-6 
and the boundary conditions 2.3-7 and 2.3-8, but the procedures can also 
be applied without modification to equations 2.3-5. The basic ideas of 
these procedures could also be applied in a solution attempt on equations 
2.3-3, however, the procedures would have to be modified to handle the 
additional relationships for r, A, A and T. 

Essentially, each of the procedures involves two solutions, an 
inner solution and an outer solution, and these must be matched. The 
outer solution satisfies the simplified disturbance equations, which are 
obtained by neglecting the terms involving mean-flow quantities since the 
mean-flow quantities are assumed to be approaching zero,and is valid for 
large values for the independent variable. The inner solution is obtained 
by integrating the disturbance equations over a range of the independent 
variable which includes most of the mean-flow boundary layer. 

The set of disturbance equations corresponding to the inner 
solution constitutes a two-point boundary-value problem for which three 
complex conditions are specified at each of the boundaries. Therefore, 
to start the integration at one end and proceed to the other end, three 
unknown complex conditions at the starting point must be guessed. Since 
the disturbance momentum equation is a fourth-order, linear, homogeneous 
differential equation, the second derivative of the disturbance stream 


function amplitude evaluated at the starting point determines the scale 
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of the disturbances and may be assigned arbitrarily. In addition to the 
remaining two unknown conditions, one complex parameter or two real para- 
meters must be determined. Thus, in total, there are six real or three 
complex eigenvalues to be determined. 

Nachtsheim [33] chose the complex wave velocity as the unknown 
parameter, and since he integrated from the surface to a suitably chosen 
edge, the other two complex eigenvalues were $¢'''(0) and S'(0). For a 
given wavenumber-Reynolds number pair, Nachtsheim obtained successive 
approximations to the eigenvalues by applying a Newton-Raphson procedure 
based on the conditions at the outer point of his integration range. To 
find a point on the neutral stability curve, he had to fix the Reynolds 
number and vary the wavenumber until the imaginary part of the complex 
wave velocity just changed sign. He then repeated this procedure at other 
Reynolds numbers to obtain the neutral stability curve. 

Knowles and Gebhart [36] followed a procedure somewhat like 
that used by Nachtsheim except they chose the real parts of the wavenumber 
and of the frequency as the unknown parameters. They set the imaginary 
part of the wavenumber to zero, as did Nachtsheim, and they also set the 
imaginary part of the frequency to zero. Therefore, for a given Grashof 
number, a point was obtained on the neutral stability curve. This procedure 
requires much less computation to obtain the neutral stability curve than 
Nachtsheim's method. Another modification introduced by Knowles and 
Gebhart was the boundary condition for the disturbance temperature differ- 
ence as discussed in section 2.3-3. 

Dring and Gebhart [37] encountered a problem with the procedures 
used by Nachtsheim and by Knowles and Gebhart, namely that the solutions 


would become unbounded if the integration proceeded far enough from the 
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surface. This effect is independent of the accuracy of the eigenvalues 
and is related to the exponential nature of the solutions in the region 
far from the surface. Theoretically the coefficients of the terms with 
positive real parts in their exponents must vanish because these solutions 
are physically inadmissible. However, numerically these coefficients are 
very small but still non-zero.. Therefore, the terms with positive real 
parts in their exponents will eventually dominate the solution if the 
integration proceeds far enough. To avoid this problem, Dring and Gebhart 
reversed the direction of integration. They used the outer solution to 
start the integration by choosing the coefficients of the terms with posi- 
tive real parts in their exponents to be identically zero. Of the remain- 
ing three complex coefficients, one coefficient was used to determine the 
scale of the solutions and the remaining two became two of the complex 
eigenvalues. In addition, they also chose the real parts of the wave- 
number and of the frequency as eigenvalues. Successive approximations to 
the eigenvalues were obtained by a Newton-Raphson procedure applied to 

the conditions at the surface. 

Although each of the above procedures appears to be relatively 
simple to apply, in practice this is not the case. Without having at 
least one complete set of eigenvalues corresponding to one point on the 
neutral stability curve, it is very difficult to guess six real or three 
complex numbers which will be close enough to the correct values for the 
iterative process to converge. This difficulty can be appreciated by 
realizing that the search technique amounts to searching for a point in 
a six-dimensional space. This problem can be overcome by noting that the 
disturbance equations can be combined to form a sixth-order, linear, homo- 


geneous differential equation which can be solved by a linear combination 
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of six linearly independent solutions. Hieber and Gebhart [39] used this 
approach while following Dring and Gebhart's procedure by integrating from 
the outer edge to the surface; but instead of guessing the coefficients 
as Dring and Gebhart did, Hieber and Gebhart integrated the solutions 
separately. They used one of the three non-zero coefficients of the linear 
combination to scale the solution and solved for the other two coefficients by 
satisfying two of the boundary conditions at the surface. Hieber and Gebhart 
used the complex wavenumber as the remaining eigenvalue. The wavenumber 
was varied until the remaining boundary condition was satisfied within 
the specified limit. This approach simplified the solution to a search 
for one complex eigenvalue rather than for three complex eigenvalues. 

The procedure followed by Hieber and Gebhart appears to be far 
superior to the procedures used previously, and therefore, this procedure 
is used in the present work. It is assumed that the disturbance amplitudes 
can be expressed as: 

> = $1 + Coby + C33 
and s = 8, + C585 + C48. 


At the starting point of the integration, the outer solutions indicate 
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where y, is the value of y at the starting point of the integration and 
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The forms of the outer solutions given in the above expression are obtained 
in exactly the same way as those obtained by Nachtsheim and by Dring and 
Gebhart. 

In their work, Hieber and Gebhart used the wavenumber as the 
complex eigenvalue, but the present work uses the real parts of the wave- 
number and of the frequency as unknown real eigenvalues. After making 
initial guesses for d- and B and separately integrating the three solutions 


from y = eto = 0, the boundary conditions 


G60) 0 
and ¢:,(0) = 0 
are used to solve for C, and C.: Then the condition on s(0) is used to 
test the approximate values for d- and Bie If this condition is not 
satisfied within the specified accuracy, e.g. 10 °, then a Newton-Raphson 
procedure is applied to the condition. For the present work s(0) = 0 is 


used as the boundary condition, and therefore, the Newton-Raphson procedure 


is of the form: 
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The partial derivatives in this procedure are approximated by finite 
differences, and \, and B. are changed to 1.001 A, and 1.001 B, to find 
these derivatives. With the new values of Ng and B, determined from this 
process, the integrations are repeated and the iteration continues until 
s(0) = 0 is satisfied within the prescribed limits. 

3.2-2 Taylor-Goertler Roll-Vortex Disturbances 

The methods of solution of the disturbance equations for the 
wave disturbances can also be applied to the equations for the roll-vortex 
disturbances. However, there are some important differences in the 
equations which result in considerable simplification. Firstly, the real 
part of B is zero since the vortices may be considered as standing trans- 
verse waves. Secondly, if the vortices are assumed to be neutrally stable 
with position and time, the imaginary parts of § and B are also zero. 
Under these simplifications the disturbance equations reduce to real 
equations. 

Equations 2.3-11, 2.3-12 and 2.3-13 represent sets of linear, 
homogeneous ordinary differential equations of eighth order. In order to 
solve these equations, a combination of eight linearly independent solutions 
is required; but following the procedure in the previous section, four of 
the coefficients must be zero because the corresponding outer solutions 
have positive real parts in their exponents. Another coefficient is set 
to unity to scale the solutions and the remaining three coefficients are 
solved for by satisfying three of the four boundary conditions given in 
2.3-14. If the above siupit cleacione are used, the equations are real and 
the remaining boundary condition is used in finding the real part of the 


wavenumber in the same way as i, and B, are obtained for the wave 


disturbance. 
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For the present work, the surface is isothermal and, therefore, 
the function a(x) is zero. This makes it possible to reduce equations 
2.3-13 to a sixth-order set, thus simplifying the procedure even further. 
Under these conditions, the disturbance amplitude functions are expressed 
as 


and Ss = 8] + C589 + C383 
At the starting point of the integration, Ye,» the solutions take the form 


of the outer solutions as follows: 
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Since err! are the only forms of the outer solutions for O51> it is 
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remainder of the procedure closely parallels that described in section 


3.2-1 and therefore will not be repeated. 
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CHAPTER IV 
RESULTS AND DISCUSSION 


4.1 Mean Flow 
foleler Errors 

Solutions of the transformed boundary-layer equations 2.2-9 
subject to the boundary conditions 2.2-10 have been obtained for air 
(Prandtl number of 0.72) adjacent to a flat, two-dimensional, isothermal 
surface. The solutions have been generated at positions along the surface 
up to X = L, where the Rayleigh number is 10®. All of the numerical 
computations have been done in double precision on an IBM 360/67 computer 
at the University of Alberta. In the computations, the value of AT was 
chosen such that Ay = 0.05 at € = 1.0 and this step size was used in all 
solutions. As a check on the truncation error in the calculations, this 
step size was doubled with the result that the eigenvalues changed in or 
beyond the third decimal place. A value of 10 * was used for &° in most 
of the calculations, but variations from 10 + to 10 © were also used with 
the result that the eigenvalues changed in the third decimal place at the 
most. For the finite-difference approximations of the partial derivatives 
with respect to &, the step size A€ was decreased until the eigenvalues 
obtained from two step sizes agreed to at least three figures for values 
of € > 0.1. The resulting step size for the first step away from € = 0 
was quite small to minimize the error in the two-point finite-difference 
approximation which is of the order of A& per step. As the solution 
proceeded in the €-direction, the value of A€ had to be increased several 
times because of a numerical stability problem. In most cases, it was 


found that the ratio of _& had to be less than about 3. Because of this 
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restriction, the last step size used in many solutions was AE = 0.4. 
This would appear to give a considerable error since the error per step 
using a three-point finite-difference formula is equal to (AE)* multi- 
plied by the third derivative of the functions with respect a Ene 
However, the solutions reveal that the functions in question are almost 
asymptotic in the vicinity of € = 1.0, and therefore, the third derivatives 
with respect to € would be very small as would be the error in the finite- 
difference approximation. The results of these computations are presented 
in figures 2 to 7. 
4.1-2 Mean-Flow Results 

Figures 2 snd 3 present typical examples of the effect of surface 
inclination on the temperature, pressure and velocity profiles at different 
positions along the surface. Figure 2 presents the results for a position 
given by € = 0.2, whereas figure 3 considers the position § = 1.0. The 
temperature and velocity profiles are shown for 0°, 30°, 45° and 90° in 
figure 2 and for 0°, 30°, 45°, 60° and 90° in figure 3. The pressure 
profiies.are.shown,for 04, 302,245.55, 600,+75 ».105,.and 120° in both 
figures. The eigenvalues for these profiles are tabulated in table l. 

Figures 4 and 5 indicate typical examples of the variations 
with position along the surface of the temperature, pressure and velocity 
profiles for specific inclinations. Figure 4 considers the profiles for 
45° eat the positions € = 0.00625, 0.01875, 0.04375, 0.1, 0.2, 0.4, 0.6 
and 1.0. Figure 5 presents the profiles for 75° at the positions € = 0.00625, 
0.025, 0.05, 0.1, 0.2, 0.4, 0.6, and 1.0. For these profiles, the eigen- 
values are listed in table 2. 

With the exceptions of the flows for 0° and 90°, the mean flow 


is non-similar and two examples of the departures from similarity of the 
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temperature, pressure and velocity are indicated in figures 6 and 7. 
Figure 6 shows this departure for a surface inclination of 45° at the 
positions €=0.00625, 0.01875, 0.04375, 0.1, 0.4, and 1.0. Figure 7 
shows the departure for an inclination of 75° at the position §=0.00625, 
0.01875, 0.1, 0.4 and 1.0. The eigenvalues for these profiles are also 
listed in table 2. 
4.1-3 Discussion of the Results 

Figures 2 and 3 clearly show that for any position along the 
surface the mean-flow boundary-layer thickness decreases and the heat 
transfer increases markedly as the surface inclination changes from 0° 
to 30°. The boundary-layer thickness continues to decrease and the 
heat transfer continues to increase at much slower rates as the in- 
clination increases from 30° to 90° at which the boundary-layer thick- 
ness reaches a minimum and the heat transfer reaches a maximum. As is 
indicated by the order-of-magnitude analysis in appendix A, the changes 
in the profiles from 0° to 30° are due primarily to the decrease in the 
lateral buoyancy force and, therefore, to the decrease in the lateral 
pressure gradient as well. For inclinations from 30° to 90°, the changes 
in the profiles are due primarily to the increase in the longitudinal 
buoyancy force. This longitudinal buoyancy force varies from zero at0° to a 
maximum at 90°, whereas the lateral buoyancy force varies from a maximum 
at.0°. to zero at,90°. 

Figures 2b and 3b indicate the variation in the pressure profiles 
as the surface varies from 0° to 120°. The pressure changes sign as the 
inclination passes through 90°, but for inclinations on opposite sides of 
the vertical, the magnitudes of the pressure are almost identical. This 


fact by itself is evidence that the pressure has little effect on the 
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flow at least over the range of inclinations from 60° to 120°. In 

support of this conclusion, the velocity and temperature profiles for 

75° and 105° as well as for 60° and 120° were found to be very nearly 
coincident. None of these profiles are included in figure 2 since the 
departures from the 90° profiles are very small. In fact for most 
purposes, the 90° profiles could be used for any angle from 60° to 120°. 
However, figure 3 does show the velocity and temperature profiles at 

60°. The departures from 90° profiles are still quite small. In this 
case, the 90° profiles could be used for angles from 75° to 105°. If 

the solutions had been obtained for surface inclinations greater than 
120°, the effect of the lateral buoyancy force should become important 
again. Then a comparison of the profiles for inclinations on opposite 
sides of the vertical should have revealed a thinner boundary layer on 
the downward-facing side, that is for a>90°. However, difficulties were 
encountered in starting the numerical procedure and the solution attempts 
had to be abandoned. The difficulties were attributed to the increasing 
importance of the lateral buoyancy force, which indirectly causes a 
positive longitudinal pressure gradient which opposes a flow in the 
positive §-direction. If the longitudinal buoyancy force is not of 
sufficient strength to overcome this positive pressure gradient, then 

it is impossible for a boundary layer to develop. It has been known 

for some time that a boundary-layer flow on a downward-facing, horizontal, 
isothermally-heated surface (a = 180°) is not possible [18, 19]. How- 
ever, it is not known if the theory will predict the break-down of the 
boundary-layer formulation at some inclination less than 180°. The 
difficulties encountered in attempting to generate solutions for inclina- 


tions greater than 120° would seem to indicate that such a limiting 
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inclination may exist, but unfortunately the numerical procedure appears 
to be incapable of finding this limit. 

Figures 4 and 5 indicate the development of the temperature, 
pressure and velocity profiles at various positions along the surface 
for surface inclinations of 45° and 75°, respectively. It is evident 
from these profiles that the thermal and momentum boundary layers con- 
tinue to grow with distance from the leading edge. This is in contra- 
diction to an assumption used in obtaining the disturbance equations 
2.3-6 and 2.3-13. However, since the neglected terms are all of the 


ea 2 
order of Ra, 5 (1+w") 


or less relative to the largest terms retained, 
this contradiction may have little bearing on the stability results. 

It has been stated previously that the boundary-layer problem 
posed for the present analysis is non-similar except for the surface 
inclinations of 0° and 90°. The departure from similarity for 45° and 
75° are indicated in figures 6 and 7, respectively. Although the magni- 
tudes of the departures indicated are a reflection of the particular 
function H(w) used in the order-of-magnitude analysis and the particular 
transformation applied to the equations, these figures do indicate 
qualitatively the departures that would result for any other choice of 
H(w), d(w) and h(w) indicated in appendices A, B and C, respectively. 

The profiles for & = 0.00625 do not appear to fit the family of curves 
indicated by the profiles for the other positions. This may be attributed, 
at least in part, to the non-zero value of & used in the lateral pressure 
gradient term to start the integration procedure as mentioned in section 
3.1. However, as mentioned in section 4.1-1, the numerical experiments 


with &° indicated that this effect was not significant for most of the 


other values of &. 
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4.1-4 Comparison with Previous Theoretical and Experimental Work 


The available information for free-convection flows over 
inclined-surfaces is still somewhat sparse. The earliest work on the 
inclined-surface problem appears to be that of Rich [22]. Rich's experi- 
mental work covered a range of angles from the vertical to 40° from the 
vertical with the heated surface facing up. For inclinations further 
from the vertical, the flow became three-dimensional. Hassan and 
Mohammed [26] attempted to explain the three-dimensional effect observed 
by Rich as a separation of the flow. However, since Rich's experiments 
were for large Grashof numbers and since he observed fluctuations in his 
results, the phenomenon might also be explained as the formation of a set 
of longitudinal roll vortices of the form observed by Sparrow and Husar 
[44]. The experimental results obtained by Rich were later found to be 
in good agreement with the theoretical and experimental results of Kierkus 
[25] and with the experimental results of Hassan and Mohammed [26]. 

Levy [23] used integral methods to study various examples of 
free-convection flow. His predictions for the heat transfer were found 
to be in good agreement with the experimental results of Rich. Levy's 
heat transfer results for a vertical and a horizontal isothermal plate 
are indicated in figures 8 and 9. These figures compare the effects of 
inclination on the heat transfer results as predicted by the present 
analysis with the results of previous theoretical and experimental investi- 
gations for € = 0.2 and for € = 1.0. The Nusselt numbers predicted by the 
present analysis are slightly higher than Levy's second approximations for 


the horizontal. For the vertical plate, Levy's second approximations are 
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much higher than the present results. However, using a quartic, Levy's 
results for the vertical plate are very close to the present results. 

The theoetical heat transfer results of Stewartson [18] and 
Ostrach [3] are also indicated in figures 8 and 9. The present results 
are in good agreement with these results. However, since the transforma- 
tion given in appendix B is based on and can be reduced to the similarity 
transformations for 0° and 90°, the results of the present analysis for 
0° and 90° would be expected to agree with the similarity solutions 
obtained by Stewartson and Ostrach, respectively. 

Michiyoshi's heat transfer results [24] are also shown in 
figures 8 and 9. The agreement between the present results and Michi- 
yoshi's results is quite good over the range of inclinations from 30° 
to 120° except near 90° for &=0.2. This effect may be due to the curva- 
ture of Michiyoshi's plate near the leading edge. The poor agreement be- 
tween 0° and 30° is probably due to the neglect of the longitudinal pressure 
gradient in Michiyoshi's work. Through its coupling with the lateral © 
pressure gradient and the lateral buoyancy force, this term is very import- 
ant for 0° since it provides the only driving force for a flow. For 
—=1.0, the zero heat transfer for 0° in Michiyoshi's work can be explained 
in the following manner. This position is the mid-point along the 
upper surface of the plate. Since the plate is finite, the flow starts 
at both ends of the plate and proceeds towards the mid-point of the 
plate, where the flow rises. in a plume. Because of the formation of a 
plume, there is no flow parallel to the surface at €=1.0. Therefore, 
from Michiyoshi's equation (6), it can be seen that the local heat 
transfer is zero. However, if the surface is semi-infinite as in the 


present analysis, then the flow starts at the leading edge and proceeds 
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along the surface without the formation of a plume. 

The heat transfer results obtained by Kierkus [25] using a 
perturbation analysis are also given in figures 8 and 9. Over the 
range of inclinations for which a comparison is possible, the results 
of the present analysis are higher than Kierkus' results, but the 
results are still in good agreement over a range of inclinations within 
about 15° of the vertical. However, the difference between the results 
increases markedly as the inclination goes beyond 15° to 20° from the 
vertical, especially for €=0.2. Since a perturbation analysis is 
limited by the size of the perturbation parameter, the increasing 
difference between the results can be explained at least partially by 
noting that Kierkus' analysis is approaching its limit of validity. 

A further comparison between the present results and Kierkus' theoretical 
results is given in table 3, in which the surface pressures are com- 
pared and found to be in good agreement. 

An interesting point about Kierkus' results is that for the 
vertical they do not agree with the vertical solutions used as his 
zeroth-order approximation (although Kierkus apparently intended this 
to be the case). Although an examination of Kierkus'equation (20) 
reveals that it is clearly non-homogeneous except for a vertical plate, 
Kierkus stated that the equation is homogeneous; on this basis, he 
concluded that homogeneous boundary conditions were unacceptable since 
they would yield only trivial solutions for the first-order approxima- 
tions. Kierkus thus established a set of non-homogeneous boundary 
conditions following a procedure used by Yang and Jerger [52], who 
analysed a free-convection flow over a vertical plate using a pertur- 


1 
bation analysis with Gr * as their perturbation parameter. In their 
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analysis, Yang and Jerger neglected terms of the order of Gr or less 
(that is second- and higher-order perturbations). Since there were no 
terms in their equations having a coefficient of Gra! all of the terms 
remaining in the equations were included in the zeroth-order approximations. 
Since the zeroth-order approximations satisfied all of the imposed boundary 
conditions of the problem, the zeroth-order solutions should have been 
the solutions to the problem at least to the order of the terms neglected. 
To improve the solutions, some of the neglected terms should be added; 
but the first-order approximations do not include any of these terms and 
therefore trivial solutions should be expected. Since Yang and Jerger 
imposed non-homogeneous boundary conditions for their first-order approx- 
imation, one is lead to believe that the boundary conditions of the 
original problem were somehow incorrect. If this were true, the zeroth- 
order approximation could have accounted for a modified boundary condition. 
It appears that the solutions presented by Yang and Jerger and by Kierkus 
are incorrect. 

Yang and Jerger stated that the boundary-layer solutions must 
match the potential-flow solutions and that one of the boundary conditions 
for their potential flow is given by their equation (21). In addition, 


(1) 


must tend to zero as y tends to infinity. 


(1) 


they stated* that y() and V 


Since the potential-flow equations are first order in y for both U and 


vy) | then yo) and yo) can each satisfy either a condition at y = 0 or 
a condition for y + © but not both. If y{) = 0 as y > © is imposed, then 
it would be fortuitous if y (1) at y = O matched the boundary-layer solution. 


Yang and Jerger, and later Kierkus, followed this approach and found that 


* Their notation is used here. 
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yl) did not match the boundary-layer solution. They modified the first- 
order boundary-layer solution accordingly, but they also imply that y6l) 
matches the zeroth-order boundary-layer solution before modification and 
that vl) = O as y >, This amounts to imposing more conditions on y (1) 
than the equations can satisfy. Therefore, it is felt that Yang and Jerger 
and Kierkus were in error in assuming that the boundary-layer solution 
follows from the potential-flow solution rather than vice versa, and 
consequently their results must be in error. 

In addition to his theoretical work, Kierkus performed some 
experiments on the free-convection flow about an inclined isothermal plate. 
He used a Mach-Zehnder interferometer to measure the velocity and temper- 
ature profiles. Hassan and Mohamed [26] also conducted an extensive 
experimental investigation of the heat transfer from an inclined, isothermal 
flat surface for a wide range of inclinations which corresponds to 0° to 
180°. They used Boelter-Schmidt heat flux meters to determine the local 
heat-transfer coefficient along the surface, and they found that most of 


their heat transfer data could be correlated within +10% by the relation 


Nu, = 0.348 (Gr, sina)}/*, 
Using this relation, some experimental points are indicated in figures 8, 
9, 10 and 11. These figures reveal that the present theoretical results 
are well within ten per cent of the above relation, thus indicating good 
agreement between theory and experiment. The theoretical results of 
Kierkus appear to be in better agreement with the experimental results, 
but in view of the above discussion of Kierkus' work, this conclusion 
may be unjustified. However, figures 8 and 9 do indicate that the present 


results are in better agreement with Hassan and Mohamed's experimental 
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Senne than are Michiyoshi's results. 

Figures 10 and 11 indicate the variation with position along 
the surface of the heat transfer for 45° and 75°, respectively. Some 
experimental results of Hassan and Mohamed are also shown. The agreement 
between theory and experiment is again very good. Kierkus' theoretical 
results are also indicated in figures 10 and 11, but again the apparent 
good agreement of his results with the present results and with those of 


Hassan and Mohamed may be fortuitous. 
4.2 Disturbance Flow 
4.2-1 Errors and Results 


Solutions of the disturbance equations have also been obtained 
for air (Prandtl number of 0.72) adjacent to a flat, two-dimensional, 
isothermal surface. For the Tollmien-Schlichting wave disturbances, 
equations 2.3-6 subject to the boundary conditions 2.3-7 and 2.3-8 have 
been solved under the conditions specified above. For the Taylor- 
Goertler roll-vortex disturbances, solutions have been obtained for 
equations 2.3-13 with the boundary conditions 2.3-14. All of the solutions 
were generated on an IBM 360/67 computer at the University of Alberta. 

The computations were done in single precision because it was found that 

a change to double precision resulted in a change in the eigenvalues in OF 
beyond the sixth significant figure. The values for Ay and y, were deter- 

mined by the values used in generating the mean-flow solutions. To obtain 
solutions of the disturbance equations at Rayleigh numbers for which there 
were no mean-flow results, the mean-flow profiles were interpolated using 

a computer-supplied routine which uses Chebyshev polynomials. For input 


parameters, this routine requires bounds on the relative loss of significance 


issokenombb-ows tel B ot insostee as 0 pee odin - 
pooaedausatb oWaw gritdotido®-noiml fol ald, 10% cooetiue | 
ot ie | 


ever 8-t.S bas T-€.S anots bm visbavod orld of “dostsive 
—rolyaT Ba! Soe | cbrede botttooga enotstbaoe orld sob | 


ror bentsado ased eva anol suloe asoied 


anortuloe sit to, aad ma S Sera bre Pi 


96 


and the relative error in the interpolation. These bounds were specified 
as 10°!° and 1078, respectively. 

Figure 12 presents the neutral stability curves for the Tollmien- 
Schlichting wave disturbances for surface inclinations of 30°, 37.5°, 45°, 
60°, 75°, 90°, 105° and 120°. Most of the curves have been generated for 
Rayleigh numbers up to 70° but the 90° curves have been extended to 
Rayleigh numbers of 2.8 x 10° on the upper branch and 4.0 x 10° on the 
lower branch. Figure 12a présents the neutral stability curves in terms 
of the wavenumber versus the Rayleigh number, and figure 12b presents the 
neutral stability curves in terms of the frequency versus the Rayleigh 
number. The curves are not extended to inclinations less than 30° because 
the boundary-layer assumptions are no longer valid in the neighborhood of 
the critical Rayleigh number. The curves do not extend beyond 120° for 
reasons noted in section 4.1. 

Figure 13 presents the neutral stability curves for the Taylor- 
Goertler roll-vortex disturbances for surface inclinations of 30°, 37.5°, 
45°, 60°, 75° and 85°. These curves indicate the variation of the wave- 
number with Rayleigh number for each surface inclination. The results are 
presented for Rayleigh numbers up to 10° or as far as convergence could be 
obtained. Inclinations at or beyond 90° are not considered because this 
mode of instability is not possible. Some of the eigenvalues obtained for 


each form of instability are tabulated in tables 4 and 5. 


4.2-2 Discussion of the Results 

The neutral stability curves for the Tollmien-Schlichting wave 
disturbances are indicated in figures 12a and 12b, and these figures are 
interpreted in the following manner. For small Rayleigh-number flows, 


there are no wavenumber-frequency combinations corresponding to a boundary- 
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layer mean flow for which the wave disturbances will grow with time or 
position. As X increases, a position is reached for which there is a 
single wavenumber-frequency pair that produces a neutrally stable distur- 
bance. For Rayleigh numbers above this "critical" Rayleigh number, there 
are ranges of wavenumbers and frequencies for which the resulting distur- 
bances grow with position or time or both. These ranges are bounded by 
the upper and lower branches of the curves shown in figure 12. The upper 
branches of the curves in figure 12a correspond to the upper branches of 
the curves in figure 12b. For any inclination, any disturbances having 
wavenumber-frequency pairs which lie above the upper branches or below 
the lower branches of bub curves in figures 12a and 12b are stable since 
they will attenuate with position or time or both. 

Figures 12a and 12b reveal that the surface inclination has a 
very significant effect on the neutral stability curves for the Tollmien- 
Schlichting wave disturbances. Since the solution for 90° was previously 
determined by Nachtsheim [33], the 90°-curves will be used as a base from 
which any changes are observed. The 90°-curves are characterized by the 
atiduence of two 'noses" and these will be referred to as the upper and 
lower noses. Nachtsheim's work revealed that in the region of the upper 
nose instability in the flow is associated with hydrodynamic effects in 
which the Reynolds stresses play a dominant role in transferring energy 
from the mean flow to the disturbances. In the region of the lower nose, 
the Reynolds stresses were found to be subtracting energy from the distur- 
bances, and the instability resulted from an energy transfer to the 
disturbances by the buoyancy effects which give rise to the temperature 
disturbances. Hence, this latter instability is characteristic of a 


thermal instability. From 90° to 30°, the lower nose becomes more 
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prominent and the region associated with the upper nose is almost completely 
absorbed. This implies that the thermal instability definitely assumes 

the dominant role as the surface is inclined further from the vertical. 
However, from 90° to 120°, the buoyancy effect is in the opposite direction 
to that for angles between 0° and 90°, and therefore it provides a stabi- 
lizing effect on the flow. 

The initial decrease in the maximum wavenumber and the maximum 
frequency as the inclination varies from 90° to 75° may be attributed to 
the decreasing role of the hydrodynamic effects noted above. Since the 
energy transfer associated with the Reynolds stresses is expected to 
continue to decrease with increasing inclination from the vertical, the 
subsequent increases in these maximum values between 75° and 30° are 
presumably due to the increasing importance of the buoyancy effects. In 
figure 12a, there is no obvious explanation for the peaks in the upper 
portion of the curves for 37.5° and 30°. Since this peak is slightly more 
prominant at 37.5°, it might be attributed to a small energy transfer 
associated with the Reynolds stresses, but the position of the peak does 
not appear to fit the pattern established for the other inclinations. 

In conjunction with the increasing effect of thermal instability 
as the surface inclination varies from 90° to 30°, there is a very signifi- 
cant increase in the range of wavenumbers associated with unstable flows. 
This wide range of wavenumbers exists even in the region close to the 
critical Rayleigh number. However, for inclinations of 45° or less, the 
ranges of frequencies associated with unstable flows decrease in the region 
close to the critical Rayleigh number. 

Figures 12a and 12b also reveal that the critical wavenumber 


increases as the inclination varies from 90° to 37.5°, but it decreases 
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slightly between 37.5° and 30° as well as from 90° to 120°. It is not 
apparent from the present analysis why, or if, the critical wavenumber 
should reach a maximum in the vicinity of 37.5°. In contrast to the 
critical wavenumber, the critical frequency is very close to a maximum 
at 90°. It decreases at 105° and 120°, but it is nearly constant from 
90° to 60° with it eas ae markedly from 60° to 30°. The critical 
frequency at 75° appears to be slightly higher than the value at 90° but 
it is difficult to determine whether or not this is due to small numeri- 
cal errors. 

The neutral stability curves for 30°, 37.5°, 45° and 60° 
indicate possible closure of the curves as the Rayleigh number tends to 
infinity. However, the profiles have not been extended far enough to 
draw any definite conclusion on this point. There is a possibility that 
the wavenumbers associated with the upper branch may tend to oscillate 
as the Rayleigh number becomes very large. This oscillatory behaviour 
of the upper branch was observed by Hieber and Gebhart [40] in their 
stability analysis of a free-convection boundary-layer flow over a vertical 
uniform-heat-flux plate. Hieber and Gebhart were unable to explain this 
phenomenon and there is nothing in the present analysis to indicate why 
such a behaviour might be observed. 

Figure 13 presents the neutral stability curves for the Taylor- 
Goertler roll-vortex disturbances in terms of the wavenumber versus the 
Rayleigh number based on the position along the surface. Unlike the 
neutral stability curves for the wave disturbances, the curves in figure 
13 do not exhibit any significant change in shape with a change in the 
surface inclination. These curves again separate the stable and unstable 


regions, but these regions are much more simply defined than for the wave 
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disturbances. For a given Rayleigh number and a particular inclination, 
any set of vortices having a wavenumber above the neutral curve will be 
damped out with time; but this set of vortices would be amplified if its 
wavenumber was below the curves in figure 13 or if the Rayleigh number is 
increased. The lower portion of each of the curves has been obtained by 
extrapolation using a five-point Adam's scheme because problems were 
encountered with the convergences of the numerical procedure in this 
region. 

The curves in figure 13 clearly show that the effect of the 
surface inclination is much greater as the inclination approaches the 
vertical. This can be seen by noting the variation in the wavenumber 
as the inclination changes for a fixed Rag eteh number. Another inter- 
esting feature of the curves in figure 13 is that each curve appears to 
be tending toward its own upper bound on the wavenumber. This reveals 
that for each surface inclination there is some size of the roll vortices 
below which the disturbances will not grow. This size increases as the 
surface inclination tends toward 90°. 

Another point which requires some consideration is the meaning 
of the critical Rayleigh number associated with the roll-vortex distur- 
bances. Following the definition for the wave disturbances, the critical 
Rayleigh number is defined as the lowest Rayleigh number at which amplifi- 
cation forat least one wavenumber just begins. For each of the curves 
shown in figure 13, the critical Rayleigh number is found by the inter- 
section of the extrapolated neutral stability curve with the zero 
wavenumber line. This implies that the roll vortices associated with the 
critical Rayleigh number are infinitely large. This is physically 


impossible, and it implies that the critical Rayleigh numbers obtained 
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in figure 13 are actually lower bounds. 

There is one other point regarding the solution of the 
disturbance equations which requires some consideration. In the 
solutions obtained for both disturbance forms, it was assumed that 
the mean flow could be approximated as a parallel flow. However, in 
section 4.1-3 it was noted that the flow is in fact non-parallel. 

It should be noted that some of the terms neglected by the parallel- 
flow assumption are of the same order as some of the terms retained 

in the equations. However, those terms retained involve the highest- 
order derivatives and the neglect of these terms would completely 
change the nature of the disturbance equations. The non-parallel 
effects are related to the inertia effects in the disturbance equations 
and if the inertia effects are considered by themselves then the terms 


at 2 
neglected are of the order of Ra 5 ELS) 


X , or less, relative to the 


terms retained. Furthermore, since the non-parallel effects are not 
associated with the highest-order derivatives, it is reasonable to 
assume that these terms can be neglected without a significant change 
in any important charateristics of the disturbance equations. 
4.2-3 Comparison with Previous Theoretical and Experimental Work 

The critical Rayleigh numbers predicted by each of the dis- 
turbance theories used in the present analysis are plotted versus the 
surface inclination in figure 14. Also shown is the theoretical critical 
Rayleigh number for a vertical, isothermal plate obtained by Nachtsheim. 
The experimental results of Lloyd and Sparrow [43] are also indicated. 
Nachtsheim's result [33] provides the only theoretical check on the 
present results for the wave disturbances. Nachtsheim used a different 


ordering analysis to obtain the set of disturbance equations for his 
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stability analysis, but the equations used in the present analysis 
have all the terms involved in Nachtsheim's equations, only the 
coefficients have changed. Since the sets of equations would be 
identical if the variables had been normalized by the same procedure, 
then the results obtained by both sets must agree within the nu- 
merical accuracy when presented in the same coordinate system. Figure 
14 reveals that this agreement has been achieved. The neutral sta- 
bility curve for 90° in figure 12a also agrees with the neutral sta- 
bility curve given by Nachtsheim. This agreement indicates that the 
numerical procedure for the wave disturbances is capable of giving 
sufficiently accurate results at least for 90°. Since this procedure 
is unaltered for other surface inclinations, there is no reason to 
doubt the accuracy of the numerical results for the other surface in- 
clinations. In fact, since the numerical procedure used for the roll- 
vortex disturbances is essentially the same as that for the wave dis- 
turbances, the above agreement is considered to be an indication of the 
accuracy of the results obtained for the roll-vortex disturbances as 
well. 

A comparison of the two disturbance theories reveals a very 
abrupt change in the form of the instability between 85° and 90° with 
another possible change between 15° and 20°. Since the latter change 
would occur at a Rayleigh number well below the limit of the boundary- 
layer approximations and since it has been obtained by extrapolation, 
this result has very little significance. In contrast, the former 
change is very significant since it implies a fundamental change in the 
mechanism of the instability. For surface inclinations of 90° or above, 


the present theory predicts that any instability in the flow initially 
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results from the formation of Tollmien-Schlichting wave disturbances 
which are two-dimensional. For surface inclinations of 85° or below, 
any instability is due initially to the formation of Taylor-Goertler 
roll vortices which make the flow three-dimensional. 

The change from one form of instability to the other at a 
surface inclination between 85° and 90° does not agree with the ex- 
perimental observations of Lloyd and Sparrow [43] which revealed a 
change between 73°and 76°. This difference is not surprising in 
view of the difference in the critical Rayleigh numbers predicted by 
theory and observed by experiment. However, in view of the work of 
Dring and Gebhart [36, 37] and Gebhart [53], the results observed by 
experiments are the final stages of the transition whereas the theor- 
etical results are predicting the initial step towards the transition 
to turbulence. For a free-convection boundary-layer flow over a uni- 
form-heat-flux vertical surface, Dring and Gebhart [37] showed that 
the disturbances which begin to amplify first do not have the highest 
amplification rate. Gebhart [53] found that the frequencies associated 
with the disturbances having the highest amplification rate were in 
good agreement with the experimental values observed by Eckert and 
Soehnghen [28], Szewczyk [32] and Lock, Gort. and Pond [42]. This agree- 
ment suggests that the difference between the theoretical and experi- 
mental critical Rayléich numbers can be attributed to an amplification 
process. 

All of the above experiments were performed without introducing 
controlled disturbances. However, using controlled disturbances intro- 
duced by a vibrating ribbon and using a Mach Zehnder interferometer to 


observe these disturbances, Polymeropoulos and Gebhart [35] performed 


eereenree re ola aeerer 
8 38 redo ott od vibtikdetant Yo sok ‘eno ont, eanado. ad? 
=xs sda Haty, 99798 ton asab "08: baw "26 mgowsed coiseadta 
8 belsovox dali’ SEAT ‘woxteq3 ad byoLs Sepemetosvecedo 


7 
i 


r 


nt a jor vat aaipapsalaian aldT tan" brett 

yd ped eae ‘erodmia dgtodied fsoktits oat nt gonsTetirb ois 2 
to axow sit 20 wolv’ ak sated daombrogxe. xe bavreedo | 
yd beyreedo ativuess ante f8el tanildeo aie tea as vaasdo0 


~kas & T9VvO | wold i doit savinastiaile om OL 
_ tara bewode [V€) dpdde bas gakvd!,so0tnwe illic | a 1098 
see sia avert ton ob er viliqms on sikged dotdw aan Py + 
betatoosss dstonsupsa3 ant, ged ine [eed icaussinllt vith r 
ok srew Stas noksaok tn tqms seodghst dita gotvest soouadausaah 
bre yx9d08 yd bovaoado aoulsv Letnomhsoa® ons soo onammorae b 
~9otge ettT . (Sd): brot, bae 2208 ,Aood baw ist rz . 


morsat esondusekb tathcdaosl gateu , sa 
| ji setsaorstisaay ‘ebnudes. goat aM Par tn 
| Bematveg [2E) ayeitds0 bab rete 


as 


gatouboriat tyorabw bout tag vs nator wnt ts we tk vad 


104 


experiments on the stability of a free-convection boundary-layer flow 
over a uniform-heat-flux vertical plate to determine the neutral sta- 
bility curve. The results of their investigation give strong support 
to the linear stability theory in predicting the initial instability. 
Hieber and Gebhart [39,40] found empirical correlations between the 
linear stability theory and the experimental observations of : a)the 
first noticeable oscillation in the boundary-layer flow, b)the first 
significant departure from a laminar flow. 

Although the above results seem to adequately account for the 
difference between the theoretical and experimental results for the 
wave disturbance, the difference for the roll-vortex disturbance re- 
mains to be explained. There does not appear to be any method of 
intoducing sets of roll vortices of a controlled size into the flow, 
and therefore, only "natural" roll-vortex disturbances can be observed. 
The theoretical predictions shown in figure 13 indicate that the flow 
initially becomes unstable due to the formation of very large roll 
vortices. In actual flows, it might be expected that the first roll 
vortices to form would be relatively small. Assuming that this is what 
actually happens, then the Rayleigh numbers indicated by the inter- 
sections of the experimental wavenumbers with the neutral stability 
curves in figure 13 could be several orders of magnitude higher than 
the values plotted in figure 14. However, since Lloyd and Sparrow did 
not determine the wavenumbers of the roll vortices first observed, the 


above explanation remains as a hypothesis. 


“ave Leviuem odd emimroasb 0% Badtg twat: 
Jxoqqua goowwe ov aotaay caoval: E 

Witideveat Iatitat ods —— 

sda(s.s 36 anoktavissde sash A ie laid vatite 

jest si3(d wok? seysl-ytebauad di Sb cobtbittose ofidnesi 

| “wold reatost 9 mov byerseqeb a if 

ef4 102 dnuosos \ioteapsbs 03. mass esfuaet evode eta: sha , 

eft xo? @iluaex Isaqemt Erion ‘bes inp vasieaualin 


“31 ssaadauialh “autov-liet edo 102, somes 


to bositzem yas od od ' ze2gqe Jon el a, 


ewolt sit otmt asta babaceeiied! 6 to seokatov Ekex ‘to: asse 
-bevagado jad aso asonsdausekb xettov-Llox “tanta” ao ot 
wolt od3.gada staotbrl EL sight ot stwode suai 
LLot sgaal yarstlo nobtemtes and ot ub) 5 dps Ay BomoN wi 
Llox taxi? ol godt boznaqus ed 4date abe sit: sa one 
senw at alkdd todd, gaimuaeéA 6. item clswttatex 8) fs O% or ” . 
—reanh odd. vd basaotbat axed secretin 
yatitdate Llaxtusn sd dakw exedmutovew, Ista 
asi tadgid sbuttagem to) e2abso Mexovee oe biboo, 
bib wouteqe bus byold hci speared ava sana at bos 
eda shavanedo terri, east axel ‘Blox ona to canceane 
‘sahestrogyii meen y 


105 


CHAPTER V 


CONCLUSIONS 


The agreement between the present theoretical and the previous 
theoretical and experimental mean-flow results appears to justify the 
boundary-layer formulation of the problem as well as the transformations 
used in the analysis. Although the solutions presented are restricted 
to constant-property free-convection boundary-layer flows of air over a 
flat, two-dimensional, isothermal surface, it can be seen from the trans-— 
formed boundary-layer equations 2.2-9 that the extension to other fluids, 
to more general two-dimensional surfaces and to more general surface- 
temperature distributions is straightforward. As pointed out in chapter 
2, it is also quite simple to extend the transformed boundary-layer 
equations to variable-property flows if the property variations are 
specified. The extension to axisymmetic flows is only slightly more 
difficult since it requires a reformulation of the Falkner-Skan and 
Lefevre transformations in the manner indicated in chapter 2. Therefore, 
using a formulation of the problem similar to that presented in chapter 
2 and using a numerical procedure such as that in chapter 3, it should 
be possible to obtain the solutions to a large sti gets of steady, lam- 
inar free-convection flows of the boundary-layer type or of the buoy- 
ant-jet type. 

The present solutions reveal that the surface inclination 
has a significant effect on the mean flow. However, the greatest effect 
of the inclination can be seen over the range from 0° to 60°. Within 
a range of about 30° on either side of the vertical, there is only a very 


small change in the velocity and temperature profiles. The profiles at 
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angles on opposite sides of the vertical are almost identical over 

this range. Over the same range of inclinations, the pressure profiles 
show a considerable change, including a change of sign as the inclina- 
tion passes through 90°; but it is apparent from the velocity and tem- 
perature profiles that the pressure variation over the boundary-layer 
has a negligible influence on the flow at least over this range of 
surface inclinations within 30° of the vertical. 

In contrast to the slight change in the mean profiles over 
the range of surface inclinations from 60° to 120°, the stability of 
the flow is very strongly influenced by the surface inclination over 
this range as reflected by the change of five to six orders of magni- 
tude in the critical Rayleigh number. Furthermore, the mechanism of the 
initial instability undergoes a change from the formation of a set of 
Taylor-Goertler roll-vortex disturbances to the formation of Tollmien- 
Schlichting wave disturbances as the surface inclination varies from 
below to above about 85°. 

For surface inclinations below approximately 80°, the present 
stability analysis reveals that any free-convection boundary-layer flow, 
which satisfies the conditions assumed in the present analysis, is 
potentially unstable if a set of sufficiently large roll vortices is 
superimposed on it. The critical size of these roll vortices decreases 
as the angle decreases. For surface inclinations below about 35°, the 
same flow is also potentially unstable if certain wave disturbances are 
introduced into the flow. The range of wavenumbers for unstable dis- 
turbances increases as the angle decreases. Therefore, assuming that 
the flow is in a region in which the boundary-layer approximations 


are valid, it may be concluded that for inclinations below about 80° the 


ovo Isoltnsbt s@omls ozs ‘leobizey joda 


#eltiorg siweesig oda aD ; 


dtsokirsv Ha 10, 
teva ealltsxq neem sda it Ag nado “saab. oy of saaanag we 
to votitdats edz ast ot 08 wor’ sotaenttonk pion, to 8 
yevo noktantioar onde oda vd boomou lant x 
-togem, 10 29 bina Kip’ oF vid 16 egress ond i . . 
ef3 io mtredsen sd3 . Pogrrosd 3 3uT rodme asl ae Bokdtas s i : 7 
fo J2se & to noltisuiot ‘als. nox? sanndo s aboays | | | 


-netmtlot 30 aokdedro? oda ‘os asonédyetebb x 


mori aeitay nolssalioak sositmue silt ag ont 2 


tnsesug of9 ,°08 visremixorggs wolsd vse ont a ie 


woll rayal-yaabaved noktosvaoo-99%3 Yrs peda) 


at ,eteyinns Jasasxg of2 nt boners bie 


7 eet 
4 
hs 


wt asokixov Lids sgtbl yldnshoizive to joe ae ae ‘seins we 
~ (it 
seaaei2sb eexitiaoy [fox gesdt i060: este Seagabaa “ 38 20) be ihreque 


etty ,°2& auods wo lsd enoltsntf{ont ah x07. ee i at 
o%6 eegnedwwIatbh svsy otadt4o I aldnsemy vifottaeiog a 


age 


a. 
Aa 
i. 

Py 


Yeou 102 exsdaunevaw) 20 22041 oat el ac 


: , spesozo9b ti pis; an onan 


co1ggs BS sti. Pitch! rn 
wolsd ariel tonbtome ‘mot Jadt aes ow 


wees , sso tetost 


107 


flow is not completely stable to all Tollmien-Schlichting wave dis- 
turbances and all Taylor-Goertler roll-vortex disturbances. However, 
for inclinations greater than 80°, there is at least some portion of 
the boundary-layer regime over which the flow is completely stable to 
all sizes of both forms of the assumed disturbances. 

The lack of theoretical and experimental stability analyses 
_with which a comparison can be made makes it difficult to draw any 
definite conclusions regarding the validity of the present results. 
However, the results do appear to be reasonable in terms of the pre- 
viously available results for 90°. In view of this apparent success 
in the present stability analysis in addition to the capability of 
generating the mean-flow solutions, it is now possible to extend the 
stability analysis to a much wider class of problems. Perhaps the 
first extension to the present work would be to determine the ampli- 
fication rates and attempt to correlate the results to the experimental 
results for both disturbance forms following procedures such as those 
used by Dring and Gebhart [37] and by Hieber and Gebhart [39,40]. 
Before proceeding with other classes of problems, the influence on 
the stability of the non-parallel effects of the mean flow should be 
determined even though the effect is expected to be small. To do 
this, the disturbance equations 2.3-6 and 2.3-13 would have to be re- 
placed by equations 2.3-5 and 2.3-12, respectively, but the numerical 
procedure should not require any modification for the wave disturbances. 
However, the procedure for the roll-vortex disturbances has to be 
modified to solve for four linearly independent solutions because of 
the additional coupling related to S,° Depending on the outcome of this 


investigation, either equations 2.3-5 and 2.3-12 or equations 2.3-6 
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and 2.3-13 could then be used to analyse the stability of constant- 
property free-convection boundary-layer flows of different fluids 
over more general two-dimensional surfaces having more general sur- 
face-temperature distributions. 

In order to consider the stability of variable-property or 
axisymmetric flows, the persinne disturbance equations 2.3-3 and 
2.3-10 can be considered for the wave disturbances and the roll-vortex 
disturbances, respectively. Although these equations are much more 
complicated, the procedure for solving these equations should be 
basically the same as the procedure used in the present work. Ob- 
taining the asymptotic solutions may be slightly more complicated, but 
once these are obtained, the procedure for obtaining the linearly 
independent solutions should be straightforward. For the wave dis- 
turbances, there will be three equations to solve instead of two 
because the disturbance stream function cannot be introduced, but 
again this should not present any serious difficulties. Therefore, 
within the limitations of the linear stability theory, the formulation 
of the equations and the procedure for solving the equations appears 
to be capable of analysing the stability of a large class of free- 


convection boundary-layer of buoyant jet flows. 
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Appendix A - A Normalization and an Order-of—-Magnitude Analysis 
of the Mean-Flow Equations 


As a first step in seeking any simplification to the set of 
steady, two-dimensional equations of motion 2.1-2, an attempt is made to 
normalize the dependent and independent variables by introducing a set 
of characteristic quantities. The definition of each of these character- 
istic quantities must come from either the boundary conditions, the 
equations or the physical description of the problem. Introducing the 
normalized variables and the characteristic quantities into equations 


2.1-2 yields the following set of equations: 
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With the equations in the above form, it is now possible to 
proceed with an order-of-magnitude analysis. In this analysis it is 
assumed that all nondimensional terms are normalized and that the order 
of each group of terms in equations A-l is determined by the set of 
characteristic quantities appearing in front of each group. 

Considering the conservation of mass equation, or the continuity 
equation, there are only two groups of terms. If the order of either group 
is much less than the order of the other group, it can be concluded that a 
special relationship exists between one of the velocities, the density 
and possibly the curvature. Since this relationship is not general, it 
follows that both groups of terms are of the same order. This implies 


the following: 


LPL: a 
ee (A-2) 


A second relationship between the characteristic quantities can 
be obtained by considering the energy equation. Since the fluid flow is 
thermally generated, it seems reasonable to expect the order of the 
advection terms to be equal to the order of the larger conduction terms, 
namely the lateral conduction terms. Equating the corresponding sets of 


characteristic quantities leads to the following relationship: 


wo = oer” (A~3) 


The remaining relationships necessary for determining the 
characteristic quantities can be found only by considering equations A-1 
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characteristic quantities found for the special cases, an attempt is made 
to generalize these quantities for any angle. 

Firstly, consider equations A-1 for a = Ox: For this surface 
orientation, the lateral body or buoyancy force in the momentum equations 
provides the only driving force for any subsequent motion. This lateral 
buoyancy force gives rise to a lateral pressure (departure) gradient, 
which in turn establishes a longitudinal pressure (departure) gradient. 
Since it is this longitudinal pressure gradient that establishes a longi- 
tudinal driving force, it follows that the lateral pressure gradient must 
be of the same order of importance as the lateral buoyancy force. This 


is expressed in the following equation: 


P 
—C = = 


In addition, the largest viscous terms in the longitudinal momentum 
equation are expected to be of the same order as the longitudinal pressure 
gradient, at least for fluids having moderate or high Prandtl numbers. 

(For low and moderate Prandtl-number fluids, the inertia terms are expected 


to be of the same order as the pressure gradient.) This leads to an expres-— 


sion of the form: 


P UU... 
=C se = 
a aay (A-5) 


Now a combination of equations A-2, A-3, A-4 and A-5 yields a set of 


characteristic quantities given by: 
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Secondly, consider equations A-1 for 0 = 1/2. For this surface 
orientation, the buoyancy force in the longitudinal momentum equation 
provides the driving force for any subsequent motion. For fluids having 
moderate or high Prandtl numbers, it is assumed that the largest viscous 
terms in the longitudinal momentum equation are of the same order as the 
buoyancy force. (For low- and moderate—Prandtl-number fluids, the inertia 


terms are assumed to be of the same order as the buoyancy force.) The 


resulting expression is: 


Hele _ P,80.8 . (A-7) 
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If it is assumed that equation A-4 applies for O=1/2 as well as for a=0, 


then a combination of equations A-2, A-3, A-4 and A-7 defines a set of 


characteristic quantities as follows: 
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We =o o Xe 
U_ = Vp Ra, *, (A-8) 
OX c 
Vo i= Yn R ‘ 
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and Po = OoVeo7 8 i 
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Finally, on the basis of the expressions given by A-6 and A-8, 
an attempt is made to extend the order-of-magnitude analysis to any angle. 
To accomplish this, consider one of the characteristic quantities, whose 
form varies with angle, and generalize this quantity by introducing an 


unknown function of the angle. For example, let 
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where H(0) = 2/5 and H(1/2) = 1/2. By inserting A-9 into equations A-2, 


A-3 and A-4, the following results are obtained: 
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If expression A-9 and A-10 are introduced into equations A-l, the following 


set of equations is obtained: 


[ + vJsina + ued + Tee 


dx 


H 
PcVc (yeu Hoy 
Ox, Seal " x (j+Jysina) 


(1+Aqy) (Yov , voy ae 
+ ( 3y + ay + AqYv} = 


2H 
OeVa® R YUrs. OU Ca yAquv 
be X “Ee lt itaqy) ox na + Tigaayy? 


1-3H 
1 oc hi R (+4/3M) ,d7u 
eta eee - "Xo ‘ees ees 


HT] 
The 
Q |< 
iN) 
2) 
i) 
PS 
Q 


1 ud?j ~- uJysino dao, ; , uJycoso, ad7a, 4 du dj 
é Gaiyaiat™ dx ni oF - dx z a: dx 
. Jycosada du, _ u als Jycosoda,* _ Ay Litas 
< dx ox (j+Jysina) 2 ‘d ; dx (1+Aqy) dx ‘dx 
dj Jycoscde OM oes u OM 
i AESEPLERNECE 2 Eee! ‘ a: ee a (fiivsinay ox 


2 OM Jycosode (1+Aqy) (M M 92y Jsina ov 
4. Saige 3 + dx)! ows 3) Laxay : (jtJysina) ox 


(M Jvcoso (m-+7 /3M) Aqav , (M+4 /3M) rAvdg 
STREET TE = (1+Aqy) ra (1+Aqy) dx 


+ 


ny bi i Vso 
Sci y \ ; ) aaa ‘ 


guiwolic2? si3 or enotseups éini boouboas — ‘OLA. om VA sotto 
; as 7 ey ii a 
‘ om " thomtsado af anot supe io 


4. 


, peboaoovl : tou + ontely ie an : oe 
+ Sb ~ xB Cont eeOREY y * 


ie to . 
0. = {vYpA ‘aaa? p 


he 
I ; 
\ ? ‘ 
FuRUDAY 4 BEVY HB wy. ae 
“(epAtL) — wo xo Cypatt) 
ate 7 alt 
al 
Hy ; : He-{ ; 
PV eer) . { yt mT ae e | B96 f es wg is ; 
a re ee ye of 2 Ea Ee) 1 el Te 
(ypAtl)” Cypaety © oF + 4%) xo Cpariy" Si. aie ‘e us 
a) uG i, b bneooylu 5 * 2B, Oakey ben. 
xb x6 KD xb 
« “ 
iu ' i ee Obs s0DY¥L) | apes, 
(vpAtt) “xb r 
r}, 
Me € ui ae 
Le oe uo MG » MGs. ahp200% 
x6’ Conta VLE Ey + Sen Bie xen. J ug oe 
a (GM | 19. ; 
evo Hotel vy" ; bie 
\ ha a 7 RT hn. 
‘S 6 (oatayt tt) © Tene! \¢ ane 
ay, ¥ 
i | f \ 
7 


‘ 4) >) : pava, (ue st + “Epa we i 
ab ok (vp: AYE) * oe . on Att a 
aie 


. bis 


123 


Jvsina dj , Jycosada A*yqv dq au 4 oM Aqu 


(jt+Jysina) Sie a dx? ~ (1+Aqy) ax! i Ox q Boe (1+Aqy) 


M2 OM av Jvsina ) , OM vy, (Pcve" Ray “Hy aM du 
8 ate Vi (EIN Pe pee ere ie eects SE etn (PAS — (—— 
oxhgys 5a oy i (jtJysina)> 4 Oy Ix y lox Se. } dy oy 


Aqu, , M,d7u Jsina du MAq _,du Jusina 
ad a eS ae en Ca Se) ee P, 
) [ay2 (jt+Jysina) By! (Teagy) ‘Oy Gtiysina): 


i MA*q7u pavar R Yq sino, 
Cukay ye lox? Ke] BO, (A-11) 
Oo Vo Ra 2H yu ov , Yvov Po Vo Ra 2H yAqu* 
Gg X te Magy ox 7 By? 7 UG X. Xo Iaewaa) 


3 , 
Pi. ab i eras 24 Ra. .O4 A PY 
a Shee od p 4 Mo, AQ ev 
Ny ay * lox, MriaheGy sa laaeh Cikday) OF 


Jsino, Ov = phil nedn co A?q?v 3M 4 dM, ov 3M 
~ (jtJysina) dy rr (jtJysina)*  (1+Aqy) 5 Gy * 3 one ; Gy 


2 OM, -Aqv Jvsina M + M d7u OM du, ,(1+Aqy) OM 
3 ay : Gtysina)! + UC 5 ay Ox By! : oe 


ts 2 2M) (oe a (1+Aqy)u @i 4 Jycosada. , , (1tAqy) + pi (ai 
3. dy “ox (jtJysina) “dx ~ dx (jt+Jysina) ‘dx 


‘ Sycosoda, , i + My du _ _MAqu (M+4 /3M) Aqu 


= dx 3’oy  (itAqy) (1+Aqy) 


est 


wp MB ; 7 ry sf e pb eee 


‘CapAeEy x6 E- oe (ypAtt) — 


i. n 
: Se 


omhevh G6) pAM le ene te MN .upAy 
ae ee abil = 
“O Ltt ve | ¢ 7 ver Cork eeta yy Pies a 


, 
{ 


“are a, get Savage tetas 
2 Mad i m—aohbinm f ; 4 Paar Hale * ~ — . 
(IL-4) ou"? : x0! . (yvpA+l) 


ws, 7 
a Ai HE 7 «| 5 t} fi .. To wy 
pei) are a «5. Voy 
= : = a - ; - [ ‘ A. agers deed = {-. fe = ore 
(VDAFL) 3 AD Wb 
4 
} 
vb em pA te Ne wn . Ma i 5a 
vo (ypAtt) v6 e [ x 
i. 3 4 : roe = 


MS. 4, Be MEd x ¥ pra.» 2 cabeitield + gs. My 


mer: path) ~ “Combat eT 
, ars ‘ | Be 
: | 


j : . ‘ , ie 


io, (ep ArE) gu ME at : f + r 


‘ 


= ee ene 4 (part) Ta 


b’ (unlwyttt) 


| 28 
oa 
anit 


i an 
: 7 a 


I 


ae 


124 


(M44 / 3M) Jusina, ;(1+Aqy) ¥ (M+4/3M) | +Jucosa da ___Aq du, 
(j+Jysina) (1+Aqy) *(jtJysina) dx (l1+Aqy) ox 


x Aqu MA udq _ 4a qdu  __Aqyu_ dq 
(1+Aqy) F isanvoee dx ox (1+Aqy) axle 
cYe R 1 OM dv , Md*v M dj 
‘ ox. =X 1 (ieaay) Ox Ox ox?  (JtJysina) Ped 


J MAy ; ¢ Y 
yeosada, > _ dq avy m (e v Ray jot cosa 


I+ 


dx (1+Aqy) dx ox Cc BO. 
k 0 Ra, Yeu (1l+a)d8 , Oda yev(1+a)d0 
! X 2 nen ox ‘ ak iu ay? 
2H 
kcO¢ OsR op 4 ven 
+f x. “Ke 5) 3 Ciraay) any) 7 
= pke%,—1 K(1+a)978 | 2K90 da , KOd*a K dj 
x. \itAay) t Ox” ve dx i ee i [Geaysina) ee 


Jycosada., - OK _ __KAy dq, [ (1+a) 90 + Oday} 
dx ox  (1+Aqy) dx Bx ds 


H 2 
kOe Ra (1+a) -Kd“8 4 OK 06 , KJsina 00 KAg 00 
tab xX ] { t39 FY a (j+Jysina) Oy x (1+Aqy) By! 
Cc 
k,.0, OsR (M+4 /3M) Ou, 5 Cun Jycosada, ” 
+ ye 1 ee aye) ox. * aoa @ : ee 
c 


yo 
Beet, a 
a 


in, 


125 


Ra ye ts i eee a + pl tar: 


SC esina) ox ‘d tS + ey 


ny J*v*sin2a j 2 (M-2/3M) | (j4Jysino.)Aqv ov , Jsinadv 
(j+Jysina)7 (j+Jysina) (1+Aqy) oy dy 


+ JAqv "sind, (M+4 /3M) 2Aqv. ou , _2Juvsina qs Jycosada 


(1+Aqy) (1+Aqy) aay) ox  (jtJysina) 2 ‘d a2 
4 —2(n2/3M) | (j4Jysina)du dv , Jvsinadu ate , Jycosada 
* “Gtdysina) Ox oy. Bed ae 


Aquv 2 Wyeosode (1+Aqy ) 2M 2Aqu 
‘ Ciskayy Ge dx) )/ i enna) “aya (+Aqy)!? 


kOe OsR M(du)* _ _2MAqu du, MA*q*ur 
+ Ux? 4 ay Gy? (1+Aqy) dy g (1+Aqy) j 


0. OsRa M dv. 2 
+ ae Xo ] (shay ox 
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The function H(wW) has no restrictions on its form as yet, though 
it seems reasonable to expect H(W) to be bounded by H(0) and H(1/2), that 
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or less relative to the largest terms in each equation are negligible, 


then equations A-1l can be simplified to the following: 
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form but it is expected to be a monotone function which must satisfy 

2/5 < H(wW) < 1/2. However, in order to obtain a numerical solution, a 
specific form must be provided for H(w). For this purpose, a reasonably 
simple, but arbitrarily chosen form is given by: 


H(w) = 2 (1+w*), (A-13) 
5 


and this is used in the present work wherever it is required. 

In order to complete the normalization and the order-of-magnitude 
analysis, the remaining characteristic quantities must be specified. For 
the properties 0, uU, and k, the properties of the fluid far from the 


surface might be used as the characteristic quantities, that is: 
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The characteristic length, Xs is considered to be the position along the 
surface at which the flow is being examined. This position is sometimes 
designated as L. The characteristic temperature difference, Ons depends 
on the surface temperature. For example, if the surface temperature is 
uniform, then 0. = oO : If the surface temperature is a known function of 
X, then the largest (or the mean) temperature difference between the sur- 
face and the undisturbed fluid might be used for 00° If there is a 


uniform heat flux along the surface, then a. is defined in terms of the 
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heat flux, i.e. 0, = GyL where q is the heat flux. Having established 
k 
these remaining characteristic quantities, the order-of-magnitude analysis 


is complete. 
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Appendix B - The Falkner-Skan Transformation 


The set of equations 2.1-3 is not easily solved even though the 
order-of-magnitude analysis has already eliminated several terms. In some 
special cases, it is possible to apply a similarity transformation and 
obtain a set of ordinary differential equations. The procedure followed 
in obtaining the similarity transformation can be extended to include a 
very wide class of problems simply by making the new dependent variables 
Enc vong of both new independent variables. This is the transformation 
referred to in Chapter I as the Falkner-Skan transformation and it is 
applicable to variable-property flows as well as to constant-property flows. 

The present work is primarily concerned with constant-—property 
flows and therefore the transformation will be derived on this basis. 
However, the extension of the procedure for variable-property flows can 
be easily seen. Starting with the set of equations 2.2-1 and assuming the 
continuity equation is eliminated by introducing a stream function, the 


following transformation is applied: 
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the angle ao The limiting values of these functions can be determined by 
considering the transformation for a = 0 and a = 1/2. When the transfor- 
mation is applied to the equations, each term in the equations contains a 
power of x. Following a procedure used in similarity transformations, the 
exponents of x in terms not involving derivatives with respect to & are 
equated. If the resulting relations lead to more than one definition of 
the exponents, then each possible transformation should be considered. If 
there are no definite advantages of one transformation over the other, 
then the choice becomes arbitrary. After choosing transformations for 
a = 0 and a = II/2, then an attempt is made to establish a general trans- 
formation for all angles. 

Consider the exponents for a = 0. One relation between the 
exponents is obtained from the conduction and advection terms in the energy 


equation, and it states that: 


RE Lee WY tee TEL EPIL S804 (B=2) 


Another relation, obtained from the pressure and buoyancy terms in the 
lateral momentum equation is: 


N +N, HF N ° : (B-3) 


A third relation, obtained from the pressure and viscous terms in the 


longitudinal momentum equation, is: 
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Solving for the other exponents gives: 


N, = nt+3 
1 5° 
N, = 4n+2 (B-6) 
5 
and N, = n-2 
4 5 


Consider the exponents for a = JI/2. It is assumed that the 
relation B-2 is also valid for a = 1/2. Then using the viscous and 
buoyancy terms in the longitudinal momentum equation, a second relation 


is: 


N, + 3N, = N,.- (B-7) 


It is also assumed that relations B-4 and B-5 hold for a = JI/2. The 


resulting set of exponents is: 


N, = nt3, 
Heat: 
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and N, = n-l 
S148C; 


A generalization of the exponents given in B-6 and B-8 is given 
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N, = nt3 , 
+ a) 
N, = 1+ n+3 + 3[nt3-d(w)] (B-9) 
d (w) d(w) 
and N, = n+3-d(w), 
d (w) 


where d(w) is any continuous function of w satisfying 
a0 JOH 55 
d(1/2) = 4 


and that 4 < d(w) < 5. 


Using the exponents given by relations B-5 and B-9 in the trans- 
formation given by B-1l and introducing the transformation into equations 


2.2-1, the following set of equations is obtained: 
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For the present work, the function d(W) has been chosen as 
dw) = 5 - 4w?. 


In the transformation B-l, it has been assumed that a stream 
function can be introduced. For some free-convection flows, e.g. certain 
axisymmetric flows, it may be impossible or inconvenient to introduce a 


stream function, in which case the transformation B-1 might be replaced 


by: 
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Appendix C - The Lefevre Tvanstocnacion 


The set of equations 2.2-6 is applicable to fluids having Prandtl 
numbers of unit order or higher. An attempt will now be made to transform 
the equations such that the new equations are applicable for all Prandtl 
numbers, including zero and infinity. . 


Consider the following transformation: 


n= 2, (0), 
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G3, 8) £,(0)m(c, E) 


and O(n, E) a MCG; Ede 


where 2,0), E, (0) and 2, (9) are disposable functions of the Prandtl 
number. Introducing this transformation into the set of equations 2.2-6 


leads to the following set of equations: 
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For the energy equation in C-1, since the flow is thermally 
generated, it is assumed that the advection terms and one of the conduction 
terms will be of the same order for the limits of 0 = 0 and 0 > ~, implying 
that CF td had Be 

aa (C-2) 
In the lateral momentum equation, it is also assumed that the pressure 
term and the buoyancy term will be of the same order for the limits of 


Oo = 0 and 0 > ©, implying that 


ma ie (C-3) 
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A third relationship can be obtained from the longitudinal momentum equation, 
but this relationship depends on a. For a = 0, it is assumed that for 
small Prandtl numbers the order of the pressure term is the same as that 


of the inertia terms, such that 
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whereas a combination of C-2, C-3 and C-5 implies that 
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A generalization of C-6 and C-7 is given by 
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terms have the same order for small Prandtl numbers, implying that 
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and combining C-2 and C-1]l leads to 
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and that 
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Substituting C-15 into C-1 leads to the following set of equations: 
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For the present work, the function h(w) has been chosen as 
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This procedure can also be applied to variable-property flows. 
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Appendix D - A Normalization and an Order-of-Magnitude Analysis 
of the Disturbance Equations 


Following the procedure outlined in appendix A for the mean- 
elo equations, an attempt is made to normalize the dependent and indepen- 
dent variables in the set of equations 2.3-1 by introducing a set of 
characteristic quantities. The definitions of these quantities must follow 
either from the boundary conditions, from the equations or from the physical 
description of the problem. Since the disturbance quantities are assumed 


to be small, it follows that: 


ope x. 
Y* << X 
Zk, <2 6X 
U * << U 
V«* << U 
W* << U 
QO * << 0 


0% “<<a 


uw * << 


~ 


me eg So) hie Vere 
cary q ay has My. Aree 
oh Bp an A 1 ; 2 rh @ io f : : 
euieas gbug 1 Lage tons gbx0 sts bus molds o xh agg 


| vie hi ces “apned aie id ‘oa: 
i¥ { Nite) sn - “a se 
as 


cise odd 26k A xtbheqas at bankiauo: lee oat gniwollo’ ’ 
-rogsbat Bas Snebesqek att $3 tteuon! fe a Biss at: squmosaa° tte enol zsups 39 

to Joa#? a anlonboxtnt vd: f-6; ¢ enotamigh ee any mk saidalgay 4 on 

wOLLG? Jeuir astyd tngup Sead; io enotatarteb it *cainnstlp staatredaas srede 
feslaydg sid woxd 76 i rokijsups aa ne s.4 (ane kAvbaos exebauod as, ony 394: f 
bomees a1 esti tinsup oome diuyekh edz sonke saottene ods to ie 989! 


= 


my! Ema) woo? at ae 


143 
Chk <0 
Cc Cc 
k* << k 
(es Cc 
pee << op 
Cc Cc 


The analysis of the mean-flow equations reveals that 


These last relationships suggest that the characteristic lengths and 
velocities for the disturbance flow might be related to the characteris-— 


tic lengths and velocities for the mean flow as follows: 
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The characteristic time, tos is chosen in such a way that the highest 

order terms involving derivatives with respect to time are of the same 
order as the highest order terms retained in each equation. The character- 
istic time and the remaining relationships between the characteristic 
quantities for the disturbance flow and the characteristic quantities for 


the mean flow may be determined by assuming that: 
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disturbance equations is obtained as follows: 
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present work, this function is 
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TABLE 1 


EFFECT OF SURFACE INCLINATIONS ON MEAN-FLOW EIGENVALUES 
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TABLE 2 


VARIATION WITH POSITION ALONG THE SURFACE 
OF MEAN-FLOW EIGENVALUES 


a = 45° 
3 sas Xr 4 
0.00000 2.0762 -0.5543  ~-4.0101 
0.00625 PIz20735 -0.4792 -2.3551 
0.01875 2.4157 -0.5717 -1.7460 
0.04375 2.9287 0.6083 -1.4169 
0.10000 3.7886 -0.6632 -1.1420 
0.20000 4.7693 -0.7194 -0.9461 
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0.00000 0.8112 -0.4146 -11.6288 
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0.01875 1.8397 0.5270 -2.2746 
0.02500 1.8995 -0.5316 -1.9479 
0.05000 2.1095 0.5411 -1.3775 
0.10000 2.3898 0.5549 -0.9645 
0.20000 2.7104 0.5733 -0.6670 
0.40000 3.0389  -0.5936 -0.4601 
0.60000 3.2144 0.6049 -0.3713 
1.00000 3.4578 -0.6207 -0.2830 
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TABLE 3 


—& = 0.2 
P ,L?o , 
ras x 10 
PRESENT KIERKUS 
-3.4306 -3.4502 
-2.5862 -2.5834 
-1.7502 -1.7370 
-0.8985 -0.8750 
— = 1.0 
Paecc , 
“pve 10 
PRESENT KIERKUS 
—-5.2694 -5.1595 
-3.9702 -3.8633 
-2.6896 -2.5976 
-1.3530 -1.3084 
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EIGENVALUES FOR THE TOLLMIEN-SCHLICHTING WAVE DISTURBANCES 


TABLE 4 


a = 120° 
Ra Ra 3 (1+w") -§ 3 (1+w*) 
x x ay 
E r Be eet 
UPPER LOWER UPPER LOWER 
BRANCH BRANCH BRANCH BRANCH 
1.000 x 10° 5.2313 1.6999 2028.3 901.4 
4.655 x 10° 3.8569 2.9651 1515.7 1253.3 
a = 105° 
1.000 x 108 5.9366 1.4709 2303.4 845.2 
3.430 x 10° 5.4234 2.6036 1930.9 1149.8 
2.383 x 10° 3.9108 1478.2 
a= 90° 
4.000 x 10° 0.8939 eos 
2.800 x 10° 12.7814 3956.4 
1.900 x 10° 7.4294 2615.7 
1.000 x 10° 7.0256 1.5044 2437.8 832.6 
1.9720: #10° 4.8850 4.6751 1587.1 1540.2 
O0e 75° 
PeUOO m 10, fp tUa7 727) al 3279) |uezo oes 744.2 
7.290 x 10° 10.9724 2917.9 
5.120 x 10° 8.8801 1.7643 OSV 2 856.3 
1.25 x 10° 7.4319 3.7765 1970.1 1263.2 
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TABLE 4 (CONTINUED) 


a = 60° 
1.000 x 10° 8.1416 1.2275 2083.3 657.5 
5.120 x 10° 10.0022 2408.3 
2.160 x 10° 12.6150 2.3081 2689.2 881.4 
6.400 x 10° 10.4028 4.3460 2062.1 1188.3 
a= 45° 
1.000 x 10° 9.5018 1.1528 1528.7 553.3 
2.160 x 10° 15.2573 2.0586 2370.5 709.9 
9,112 x 10° 18.5805 2894.0 
8.000 x. 10° 15.9489 7.8290 1797.5 2105.1 
a = 37.5° 
1.000 x 10° 9.5015 1.0921 1393.6 484.8 
2.160 x 10° 15.0876 1.8732 2146.5 601.3 
6.400 x 10° 20.3577 2.8793 2916.4 698.5 
8.000 x 10° 20.9067 5.4852 2155.6 785.5 
1.000 x 10° 18.4608 10.1014 1301.7 813.4 
5.12 x 10* 16.8832 12.2704 1059.4 824.0 
a@ = 30° 
1.000 x 108 9.4161 1.0081 1280.0 404.2 
2.160 x 10° 15.0306 2009.3 
6.400 x 10° 20.8175 2.4040 2756.3 535.0 
4.288 x 10° 22.6427 2974.3 
8.000 x 10° 24.5873 4.1110 2419.9 557.3 
1.000 x 10° 23.7006 6.5226 1570.8 519.4 
oh 16.1803 12.5260 679.6 529.4 
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TABLE 5 


EIGENVALUES FOR THE TAYLOR-GOERTLER ROLL-VORTEX DISTURBANCES 
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Ray z Re Sor 

OlS657 4 O75" 1 Ore 605- a = 945° sate 37507 ase 1305 
1.000 x 10° 12.920 23.896 32.098 
Sys Dba al 11.970 23.301 
2.160 x 10° 100277 22.415 ~~ =+=30-960 . 36.186 *: 38.176 
6.400 x 10° 6.668 20.933 29.725 35.246 37.341 39.087 
eT 0ux tor 1.915 | 
8.000 x 10° 17.784 26.964 33.103 35.457 37.426 
1.060 x 10° 13.484 23.064 30.063 32.645 34.795 
1.250 x 107 6.956 
83.74 16.459. 24.533 227.275 29.524 
30.52 139515) 215429) 245248 | 926.379 
6.592 16.059 18.343 20.291 
1.953 102778 126425 4 “3.935 
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